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Abstract 



of freedom. The microsystem evolution is then described by tracing over the 



Ft , , Dissipation, the irreversible loss of energy and coherence, from a microsys- 

tem, is the result of coupling to a much larger macrosystem (or reservoir) 
^ \ which is so large that one has no chance of keeping track of all of its degrees 

o 

o 

<N . 

f^ , reservoir states, resulting in an irreversible decay as excitation leaks out of the 

^N \ initially excited microsystems into the outer reservoir environment. Earlier 

^, 
CIh, treatments of this dissipation described an ensemble of microsystems using 

-(— > ' 

G ' density matrices, either in Schrodinger picture with Master equations, or in 

^ , 

qh, Heisenberg picture with Langevin equations. The development of experimen- 



tal techniques to study single quantum systems (for example single trapped 



H \ ions, or cavity radiation field modes) has stimulated the construction of theo- 

retical methods to describe individual realizations conditioned on a particular 
observation record of the decay channel, in the environment. These methods, 
variously described as Quantum Jump, Monte Carlo Wavefunction and Quan- 
tum Trajectory methods are the subject of this review article. We discuss their 
derivation, apply them to a number of current problems in quantum optics 
and relate them to ensemble descriptions. 
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I. INTRODUCTION 

Quantum mechanics is usually introduced as a theory for ensembles, but the invention 
of ion traps , for example, offers the possibility to observe and manipulate single particles, 
where observability of quantum jumps, which are not be seen directly in the ensemble, lead 
to conceptual problems of how to describe single realizations of these systems. Usually Bloch 
equations or Einstein rate equations are used to describe the time evolution of ensembles of 
atoms or ions driven by light. New approaches via conditional time evolution, given say when 
no photon has been emitted, have been developed to describe single experimental realizations 
of quantum systems. This leads to a description of the system via wave functions instead 
of density matrices. This conditional "quantum trajectory" approach is still an ensemble 
description, but for a sub ensemble where we know when photons have been emitted. 

The jumps that occur in this description can be considered as due to the increase of 
our knowledge about the system which is represented by the wave-function (or the density 
operator) describing the system. In the formalism to be presented one usually imagines 
that gedanken measurements are performed in a rapid succession , for example, on the 
emitted radiation field. These will either have the result that a photon has been found in 
the environment or that no photon has been found. A sudden change in our information 
about the radiation field (for example through the detection of a photon emitted by the 
system into the environment) leads to a sudden change of the wave-function of the system. 
However, not only the detection of a photon leads to an increase of information but also 
the failure to detect a photon (i.e. a null result). New insights have been obtained into 
atomic dynamics and into dissipative processes, and new powerful theoretical approaches 
developed. Apart from the new insights into physics these methods also allow the simulation 
of complicated problems, e.g., in laser cooling that were completely intractable using the 
master equation approach. In general it can be applied to all master equations that are of 
Lindblad form, which is in fact the most general form of a master equation. 

This article reviews the various quantum jump approaches developed over the past few 



years. We focus on the theoretical description of basic dynamics and on simple instructive 
examples rather than the application to various numerical simulation methods. 

Some of the topics covered here can also be found in earlier summaries (Cook, 1990; Erber 
et ai, 1989; M0lmer and Castin, 1996) and more recent summer school lectures (Knight and 
Garraway, 1996; M0lmer, 1994; ZoUer and Gardiner, 1995). 

II. INTERMITTENT FLUORESCENCE 

Quantum mechanics is a statistical theory which makes probabilistic predictions of the 
behaviour of ensembles (an ideally infinite number of identically prepared quantum systems) 
using density operators. This description was completely sufficient for the first 60 years of the 
existence of quantum mechanics because it was generally regarded as completely impossible 
to observe and manipulate single quantum systems. For example, Schrodinger, in 1952 wrote 
(Schrodinger, 1952). 

... we never experiment with just one electron or atom or (small) molecule. In thought- 
experiments we sometimes assume that we do; this invariably entails ridiculous consequences. 
{...} In the first place it is fair to state that we are not experimenting with single particles, 
any more than we can raise Ichthyosauria in the zoo. 

This (rather extreme) opinion was challenged by a remarkable idea of Dehmelt which he first 
made public in 1975 (Dehmelt, 1975; Dehmelt, 1982). He considered the problem of high 
precision spectroscopy, where one wants to measure the transition frequency of an optical 
transition as accurately as possible, e.g., by observing the resonance fiuorescence from that 
transition as part (say) of an optical frequency standard. However, the accuracy of such a 
measurement is fundamentally limited by the spectral width of the observed transition. The 
spectral width is due to spontaneous emission from the upper level of the transition which 
leads to a finite lifetime r of the upper level. Basic Fourier considerations then imply a 
spectral width of the scattered photons of the order of r^^. To obtain a precise value of the 
transition frequency, it would therefore be advantageous to excite a metastable transition 



which scatters only a few photons within the measurement time. On the other hand one 
then has the problem of detecting these few photons and this turns out to be practically 
impossible by direct observation. So obviously one has arrived at a major dilemma here. 
Dehmelt's proposal however suggests a solution to these problems, provided one would be 
able to observe and manipulate single ions or atoms which became possible with the invention 
of single ion traps (Paul et al, 1958; Paul, 1990) (for a review see (Horvath et al, 1997). We 
illustrate Dehmelts idea in its original simplified rate equation picture. It runs as follows. 

Instead of observing the photons emitted on the metastable two-level system directly, he 
proposed to use an optical double resonance scheme as depicted in Fig. 1. One laser drives 
the metastable <-* 2 transition while a second strong laser saturates the strong ^^ 1; the 
lifetime of the upper level 1 is for example 10~^s while that of level 2 is of the order of Is. 
If the initial state of the system is the lower state then the strong laser will start to excite 
the system to the rapidly decaying level 1, which will then lead to the emission of a photon 
after a time which is usually very short (of the order of the lifetime of level 1). This emission 
restores the system back to the lower level 0; the strong laser can start to excite the system 
again to level 1 which will emit a photon on the strong transition again. This procedure 
repeats until at some random time the laser on the weak transition manages to excite the 
system into its metastable state 2 where it remains shelved for a long time, until it jumps 
back to the ground state either by spontaneous emission or by stimulated emission due to 
the laser on the <-^ 2-transition. During the time the electron rests in the metastable 
state 2, no photons will be scattered on the strong transition and only when the electron 
jumps back to state can the fluorescence on the strong transition start again. Therefore 
from the switching on and off of the resonance fluorescence on the strong transition (which 
is easily observable) we can infer the extremely rare transitions on the ^^ 2 transition. 
Therefore we have a method to monitor rare quantum jumps (transitions) on the metastable 
^^ 2 transition by observation of the fluorescence from the strong ^^ 1 transition. A 
typical experimental fluorescence signal is depicted in Fig. 2 (Thompson, 1996) where the 
fluorescence intensity I(t) is plotted. However, this scheme only works if we observe a 



single quantum system, because if we observe a large number of systems simultaneously the 
random nature of the transitions between levels and 2 implies that some systems will be 
able to scatter photons on the strong transition while others are not because they are in their 
metastable state at that moment. From a large collection of ions observed simultaneously 
one would then obtain a more or less constant intensity of photons emitted on the strong 
transition. The calculation of this mean intensity is a straightforward task using standard 
Bloch equations. The calculation of single system properties such as the distribution of the 
lengths of the periods of strong fluorescence, required some effort which eventually led to 
the development of the quantum jump approach. Apart from the interesting theoretical 
implications for the study of individual quantum systems, Dehmelt's proposal obviously has 
important practical applications. An often cited example is the realization of a new time 
standard using a single atom in a trap. The key idea here is to use either the instantaneous 
intensity or the photon statistics of the emitted radiation on the strong transition (the 
statistics of the bright- and dark periods) to stabilise the frequency of the laser on the weak 
transition. This is possible because the photon statistics of the strong radiation depends 
on the detuning of the laser on the weak transition (Kim et ai, 1987; Kim and Knight, 
1987; Kim, 1987; Ligare, 1988; Wilser, 1991). Therefore a change in the statistics of bright 
and dark periods indicates that the frequency of the weak laser has shifted and has to be 
adjusted. However, for continuously radiating lasers this frequency shift will also depend on 
the intensity of the laser on the strong transition. Therefore in practise pulsed schemes are 
preferable for frequency standards (Arecchi et al, 1986; Bergquist et al, 1994) 
Due to the inability of experimentalists to store, manipulate and observe single quantum 
systems (ions) at the time of Dehmelt's proposal, both the practical as well as the theoretical 
implications of his proposal were not immediately investigated. It was about ten years later 
that this situation changed. At that time Cook and Kimble published a paper (Cook and 
Kimble, 1985) in which they made the first attempt to analyse the situation described above 
theoretically. Their advance was stimulated by the fact that by that time it had become 
possible to actually store single ions in an ion trap (Paul trap) (Paul et al, 1958; Neuhauser 
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et ai, 1980; Paul, 1990). 

In their simplified rate equation approach Cook and Kimble started with the rate equa- 
tions for an incoherently driven three level system as shown in Fig. 1 and assumed that 
the strong ^^ 1 transition is driven to saturation. They consequently simplify their rate 
equations introducing the probabilities P+ of being in the metastable state and P of being 
in the strongly fluorescing <-^ 1 transition. This simplification now allows the description 
of the resonance fluorescence to be reduced to that of a two state random telegraph process. 
Either the atomic population is in the levels and 1 and therefore the ion is strongly radi- 
ating (on), or the population rests in the metastable level 2 and no fluorescence is observed 
(off). They then proceed to calculate the distributions for the lengths of bright and dark 
periods and find that their distribution is Poissonian. Their analysis that we have outlined 
very briefly here is of course very much simplified in many respects. The most important 
point is certainly the fact that Cook and Kimble assume incoherent driving and therefore 
adopt a rate equation model. In a real experiment coherent radiation from lasers is used. 
The complications arising in coherent excitation finally led to the development of the quan- 
tum jump approach. Despite of these problems the analysis of Cook and Kimble showed the 
possibility of direct observation of quantum jumps in the fluorescence of single ions, a pre- 
diction that was confirmed shortly afterwards in a number of experiments (Bergquist et ai, 
1986; Nagourney et ai, 1986a; Nagourney et ai, 1986b; Sauter et ai, 1986b; Sauter et ai, 
1986c; Dehmelt, 1987) and triggered off a large number of more detailed investigations start- 
ing with early works by Javanainen in (Javanainen, 1986a; Javanainen, 1986c; Javanainen, 
1986b). The following effort of a great number of physicists eventually culminated in the 
development of the quantum jump approach. Before we present this subsequent develop- 
ment in greater detail we would like to study in slightly more detail how the dynamics of 
the system determines the statistics of bright and dark periods. Again assume a three-level 
system as shown in Fig. 1. Provided the ^^ 1 and ^^ 2 Rabi frequencies are small 
compared with the decay rates, one finds for the population in the strongly-fluorescing level 
1 as a function of time something like the behaviour shown in Fig. 3 (we derive this in detail 
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in a later section.). We choose for this figure the values 71 ^ 72, for the Einstein-coefficients 
of levels 1 and 2 refiecting the metastability of level 2. For times short compared with the 
metastable lifetime 72"^, then of course the atomic dynamics can hardly be aware of level 2 
and evolves as a - 1 two-level system with the "steady state " population pn of the upper 
level. After a time 7^^, the metastable state has an effect and the (ensemble- averaged) 
population in level 1 reduces to the appropriate i/iree- /ewe/ equilibrium values. The "hump" 
Apii shown in Fig. 3 is actually a signature of the telegraphic fluorescence discussed above. 
To show this, consider a few sequences of bright and dark periods in the telegraph signal as 
shown in Fig. 4. The total rate of emission R is proportional to the rate in a bright period 
times the fraction of the evolution made up of bright periods. This gives 

but this has to be equal to the true average, 

i? = 7iPii(oo), (2) 

so that 



Td Pii-Pii(oo) Ap 
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(3) 



Tb Pii(oo) Pii(oo) 

and the ratio of the period of bright to dark intervals is governed, as we claimed, by the 
"hump" Apn. 

So far, we have concentrated on situations where the Rabi frequencies are small (or 
for incoherent excitation). What happens for co/iereni resonant excitation with larger Rabi 
frequencies? The answer to this question is nothing (Knight et ai, 1986): there are essentially 
no quantum jumps, at least at any signiflcant level for coherently-driven resonantly excited 
three-level systems! But this is because of the idea of resonance is tricky: the strong Rabi 
frequency on the <-> 1 transition dresses the atom and the AC-Stark effect splits (Autler 
and Townes, 1955; Knight and Milonni, 1980) the transition, forcing the system substantially 
out of resonance. If this is recognised and the probe laser driving the ^^ 2 transition is 
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detuned from the bare resonance until it matches the dressed atom resonance, then the 
jumps and telegraphic fluorescence return. We investigate this phenomenon more closely in 
Section V. As far as we know, the dependence of the telegraph fluorescence on detuning for 
coherently excited transitions has yet to be confirmed experimentally. 

Let us return to the idea of a null measurement. We imagine that we observe the 
fluorescence from a driven three-level ion over a time scale which is long compared with the 
strongly fluorescing state lifetime 7]"^ but very short compared with the shelf state lifetime 
72^ , so that 7f ^ At ^ 72" . Pegg and Knight (Pegg and Knight, 1988a; Pegg and Knight, 
1988b) have shown that the average period of brightness and darkness in the telegraphic 
fluorescence can be obtained very straightforwardly from considerations of null detection. 
During such an interval At, the population in the shelf state, P2{t) hardly has time to evolve, 
but population can be rapidly cycled from the ground state |0) to the strongly-fluorescing 
state |1) and back. Detection of a photon at the beginning of a At interval implies a survival 
in the 0-1 sector for the whole interval and a bright period, whereas a null detection is 
sufficient for us to be confident that the atom is shelved for the whole At interval and a 
dark period ensues. 

If we take our origin of time to be after an interval At in which we see a photon, then 
-^2(0) = 0. We can introduce the "life expectancy" Tb as the time the atom spends in the 
0-1 sector continuously. If the atom is still in this sector at a time ti (known from an 
observation of another fiuorescence photon just prior to ti), then the life expectancy will 
be also be Tg. So we can partition the outcomes into the case where at ti it has survived 
in the 0-1 sector with probability Pio(ti) and the case where the ion did not survive the 
whole interval ti continuously in the 0-1 sector (Pegg and Knight, 1988a; Pegg and Knight, 
1988b) 

Tb = Pioih) (ti + Tb) + (1 - Pio(ti))/ti , (4) 

where / is a fraction (< 1). Then for small ti 

^^ " 1 - Pio(ti) " 1 - Pio(ti) '' ' ^^^ 
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given P2(0) = 0. (6) 

i=0 



and if ti is small so we may neglect the possibility of a return from state |2) back in to the 
0-1 sector, (1 - Pio(ti)) ^ ^^2(^1), so that 

This is finite, so we know that the fluorescence will terminate. To obtain a value for Tg, we 
merely need to know the evolution equation (not its solution) for the population in state 
1 2) : this would be the Bloch equation for coherent excitation, or the Einstein rate equation 
for incoherent excitation. 

The calculation of the mean period of darkness proceeds along similar lines: if no photons 
are detected in an interval At just before t = 0, we find 



^^ " dt 



given P2(0) = l. (7) 

t=o 



The analysis presented here obviously also applies for the density operator equations in 
exactly the same form and we obtain 

TE' = {P22\=o Withp22(0)=0 , (8) 

TB' = -{P22\=0 Withp22(0) = l . (9) 

Here the dot means the average gradient of the P22 versus t curve over a range of order 
At. Because At is much smaller than the characteristic change of P22 it is very close to the 
normal derivative at all points. 

It is straightforward to use this idea of "collapse by non-detection" to estimate the char- 
acteristic time needed to be sure that a quantum jump has occurred (Pegg and Knight, 
1988b). There are Ts/td times as many short dark periods between photon emissions as 
there are prolonged dark periods of average length To, where td (^ 7]"^ for strong transition 
saturation) is the average length of the short period. Thus the probability that an emis- 
sion will be followed by a long dark period is approximately td/Ts for Tg ^ t^^, and the 
probability that it will be followed by a short dark period is close to unity. 

Immediately following a photon emission, a dark period of length at least r (with r < Tp) 
can exist for two complementary reasons: (a) the atom goes to state |2) and therefore does 
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not decay for a time of the order of T^, or (b) the atom is still in the |0) ^^ |1) plane but has 
not yet emitted a photon. The probability of (a) occurring is t^/TB and the probability of 
(b) is approximately exp{—T/td). Clearly for r < t^^ it is much more likely that any observed 
dark period of length r is due to (b), but this becomes rapidly less likely as r increases. The 
point where the observation of the dark period is just as likely to involve (a) as (b) is found 
by equating the two expressions to give 

^ = ^ HTbIU) . (10) 

J-B J-B 

It follows that the sampling period At must be greater than r given by Eq. (10) this value 
in order for the observation of darkness during At to imply (a) with reasonable certainty. 

Further, because we know that immediately following the emission the atom is in |0), so 
the probability of being in |2) is zero, and because Eq. (10) gives the order of the time of 
darkness required for the probability of being in |2) to grow to about |, Eq.(lO) gives the 
characteristic time for the wave-function collapse by non-detection. This characteristic time 
can be associated with the time necessary for us to be certain that a quantum jump from |0) 
to |2) has occurred. For completely coherent excitation Eq. (10) reduces to an expression 
similar to that for the "shelving time" found in (Porrati and Putterman, 1987) and in (Zoller 
et a/., 1987). Note that while the collapse by non-detection of the system into the metastable 
state requires a finite time the collapse of the wave function due to the detection of a photon 
has to be viewed as practically instantaneous. When we detect a photon our knowledge of 
the system changes suddenly and this sudden change of knowledge is reflected by the sudden 
change of the system state which, after all, represents our knowledge of the system. 

III. ENSEMBLES AND SHELVING 

Before we develop detailed theoretical models to describe individual quantum trajectories 
(i.e. state evolution conditioned on a particular sequence of observed events), it is useful to 
examine how the entire ensemble evolves. This is in line with the historical development as 
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initially it was tried to find quantum jump characteristics in the ensemble behaviour of the 
system. We do this in detail for the particular three-level V-configuration (shown in Fig. 1) 
appropriate for Dehmelt's quantum jump phenomena. For simplicity, we examine the case of 
incoherent excitation. Studies for coherent excitation using Bloch equations can be found 
, for example, in (Kim, 1987; Kimble et ai, 1986; Ligare, 1988; Nienhuis, 1987; Schenzle and 
Brewer, 1986). The Einstein rate equations for the V-system are (Pegg et ai, 1986b) 

j^Pn = -iAi + BiWi)pn + BiWipoo , (11) 

-^P22 = -(^2 + B2W2)P22 + B2W2P00 , (12) 

-Poo = -{BiWi + B2W2)poo + {Ai + BiWi)pu + (^2 + 53^^2)^22 , (13) 

where A^, Bi are the Einstein A- and B-coefiicients for the relevant spontaneous and in- 
duced transitions, Wi the applied radiation field energy density at the relevant transition 
frequency, and pu is the relative population in state i (poo + Pii + P22 = 1 for this closed 
system). In shelving, we assume that both BiWi and Ai are much larger than B2W2 and 
A2, and furthermore that B2W2 ^ ^2- The steady-state solutions of these rate equations 
are straightforward to obtain, and we find 

( ^ B^W^iA2 + B2W2) 

Pill ^ooj Ai{A2 + 2B2W2) + BiWi{2A2 + 3B2W2) ' ^ ^ 

.^ ^ ^ B2W2iA, + B,W,) 

P22I ^ooj Ai{A2 + 2B2W2) + BiWi{2A2 + 3B2W2) ■ ^ ^ 

Now if the allowed ^^ 1 transition is saturated 

pn(t ^ 00) ^ , o D TT/ ~ Poo{t -^ 00) , (16) 



2A2 + 3B2W2 



and 



B2W2 1 

P22(t^oo)^ ^^^^3^^^^ ^- . (17) 

Now we see that a small B2W2 transition rate to the shelf state has a major effect on the 
dynamics. Note that if the induced rates are much larger than the spontaneous rates, the 
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steady state populations are poo = Pii = P22 = | : the populations are evenly distributed of 
course amongst the constituent states of the transition. 

However, the dynamics reveals a different story to that suggested by the steady state 
populations. Again, if the allowed transition is saturated, then the time-dependent solutions 
of the excited state rate equations tell us that for poo(O) = 1 we have 

„ (i-\ ^ B2W2 „-(A2+3B2W2/2)t _ }_ -{2BiWi+Ai+B2W2/2)t , ^2 + ^2^2 f-.^. 

^''^' 2(2^2 + 3521^2) 2 ^2A2 + 3B2W2 ' ^ ' 

and 

Note that the these expressions are good only for strong driving of the ^^ 1 transition. 
This especially means that for short times poo is of the order of 1/2 which results in Eq. (19). 
Then for a very long-lived shelf state 2, we see that for saturated transitions {BiWi ^ Ai) 

Puit) ^l{l + l (e-^^^^^*/2 - 3e-^^^^^*)} , (20) 



and 



P22 



(t)^-{l-e-3S^^^*/2} (21) 



These innocuous-looking expressions contain a lot of physics. We remember that state 1 is 
the strongly fluorescing state. On a time-scale which is short compared with (-021^2)"^, we 
see that the populations attain a quasi-steady state appropriate to the two-level (0 ^^ 1) 
dynamics 



Pi 



i(t short) ~i{l-e-2^i^i*}^i . (22) 



This can of course be conflrmed in an experiment (Finn et ai, 1986; Finn et ai, 1989). For 
truly long times the third, shelving, state makes its effect and 

pn(tlong)~^ , (23) 
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as we saw qualitatively in Fig. 3. As we saw earlier in the discussion of Eq. (3), the 
change from two-level to three-level dynamics already gives us a signature of quantum 
jumps and telegraphic fluorescence provided we are wise enough to recognise the signs. 
Figure 3 illustrates the change from two to three-level dynamics. 

The steady-state populations are sufficient to describe the average level of the fluorescent 
intensity. But how do quantum jumps and shelving show up in the intensity correlations? 
For example, let us examine the second order correlation function 

where the colons describe normal ordering (Loudon, 1983). This correlation function is 
straightforward to compute from the Einstein rate equations using the quantum regression 
theorem (Lax, 1963) which relates one-time to two-time correlations given the dynamics 
is Markovian. Here of course there are two intensities, that of the 1^0 and of the 
2^0 fluorescence, so we can correlate the two light fields: "1" with "1", or "1" with 
"2" and so on, where "1" and "2" represent the fiuorescence on the 1^0 and the 2 — i> 
transitions respectively. So let us concentrate on evaluating glj{t,T), which represents the 
joint probability of detecting a fiuorescent photon j (j = 1, 2) on the transition j at time t 
and some other photon (not necessarily the next photon) from transition i at time (t + r). 
it is straightforward to show (Pegg et ai, 1986b) that 

^ B2W2 (A2+3B^W2mT _ 2^2 + 3B2W2 (2BrWj+A^+B2W2/2)r (2^^ 

^2{A2 + B2W2) 2{A2 + B2W2) ' ^ ' 

and 

9?^ir) = g'S{r) = 1 - e-^^^^'^^'^^/'^ . (26) 

Using BiWi ^ B2W2 and that the transitions are saturated, then these correlation functions 
simplify to give 

9?^{r) = ^!?(r) = 1 + 1 (e-^^^^-/^ - 36"^^^^-) , (27) 
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so that as expected, the correlation functions obey the same evolutions as the populations. 
It is worth noting that for short times r we expect to see anti-bunching (Loudon, 1983) from 
this three— level fluorescence, and this has been observed experimentally from trapped ions 
(Itano et al, 1988; Schubert et al, 1992). 

In the mid 1980's when studies of quantum jump dynamics of laser-driven three-level 
atoms began in earnest, a great deal of effort was expended in determining the relationship 
between joint probabilities of detection of a photon at time t and the next or any photon 
a time r later. This was addressed in detail by Cohen- Tannoudji and Dalibard (Cohen- 
Tannoudji and Dalibard, 1986; Reynaud et al, 1988), by Schenzle and Brewer (Schenzle et 
al., 1986; Schenzle and Brewer, 1986) and others (Cook, 1981; Lenstra, 1982). One attractive 
approach, advocated by Cohen- Tannoudji and Dalibard, uses a dressed manifold and from 
this evaluates the delay function describing the distribution of delay times before the next 
emission occurs. In Fig. 5 that laser excitation couples together these atom+field states, 
but fluorescence does not: spontaneous emission is an irreversible loss out of this manifold 
to states with reduced photon number in the excitation modes, but with photons created in 
initially un-occupied free-space modes. 

We expand our atom+fleld state vector into the basis states |/) with flxed number of 
fluorescence photons in the radiation fleld as shown in Fig. 5 as 

|^(t)) = ^aKt)e-^'*/^|/) , (28) 

I 

and solve for the probability amplitudes ai{t). The probability of remaining without further 

emission in the n-excitation atom+fleld manifold up to time r is 

Po{r) = j:W^{T)\' . (29) 

i 

The negative differential of this survival probability describes the delay function /i(r) 
(Cohen- Tannoudji and Dalibard, 1986) 

/.(.) . -f . (30) 

15 



so that the probabihty of there being an interval r between one photon being emitted 
(detected) and the next is 

Po(r) = l- r hir')dr' . (31) 

Jo 

To evaluate this interval distribution function it is sufficient to calculate the atom+field 

dressed state amplitudes. This demonstrates the utility of this approach: there is no need to 

solve the potentially complicated Bloch equations for the driven three-level atom, although 

of course this can be done (Schenzle and Brewer, 1986). Kim and coworkers (Kim et ai, 

1987; Kim and Knight, 1987), Grochmalicki and Lewenstein (1989a), Wilser (1991) and 

later others have used the delay function to describe the shelving in the V system. The 

conditional probability of an atom emitting any photon between time r and t + dr after 

emitting a photon at time r = is I{T)dT. The photon emitted at time r can be the first 

to be emitted after that at r = or the next after any one at time r' (0 < r' < r), so that 

I{r) = h{r)+ rdr'l{r')h{r-r') , (32) 

JO 

where Ji (r) is the interval distribution. If this is expressed in terms of its Laplace transform 
/(z), then 

1-/1(2;) 

The function g^'^^r) is the normalized correlation function for a photon detection at r = 
followed by the detection of any photon (not necessarily the next) at time r. It follows that 

i r ^ 00 



or equivalently in Laplace space 



^(2)(^) = M.__ , (35) 



so that 



9^'\z) 



hm.^-^^^^) 



^0 zh{z) 
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''^'^ (36) 



l-h(z) 



For the case of incoherent, rate-equation excitation of a three-level V-system atom, it is 
straightforward to calculate the atom+field survival probability, differentiate this to generate 
the delay function and from Eq. (36) deduce g^'^^r). If this is done, precisely the same form 
is obtained as that from the quantum regression theorem. The merit of this approach is 
easier to appreciate for the case of coherent excitation, where the regression theorem requires 
the solution of the three-level Bloch equations and the solution of 8th order polynomial 
characteristic equations, compared with the need to solve three coupled equations, using the 
delay function route. In Section V we further illustrate the connection between the "next" 
photon probability density and the "any' photon rate in Eqs. (183)-(194). 

Rather than examine the correlation functions, it may be useful to examine other prop- 
erties of the photon statistics, and in particular the variance in the photon numbers in the 
detected fluorescence (Kim et al., 1987; Kim and Knight, 1987; Kim, 1987; Jayarao et al., 
1990). For a Poissonian field, (An)^ = (n), but for a sub-Poissonian field (An)^ < (n), 
and for a super-Poissonian field, (An)^ = (n). To characterise the deviation from pure 
Poissonian fluctuations, Mandel (Mandel, 1979) defined the parameter 

,,, _ (An)^ - (n) 
QM{rj = T-r , (37j 

which can be written in terms of the mean intensity (/) and the second order correlation 
function g'^'^\t) as 

QM{r) = ^ {[dt2fj dt,g^'\t,)^ - {I)t . (38) 

If this is used to describe the fluorescence from a three-level atom involving shelving (Kim 
et al, 1987; Kim and Knight, 1987; Kim, 1987) then as r ^ oo for saturated transitions 

gM(r^ oo) = ^(/) , (39) 

-LB 

where Tb,Td are the mean times of bright and dark periods in the telegraphic fluorescence 
signal from the three-level atom shown in Fig. 1. If a dark period does not occur, Qm{t — >■ 
oo) -^ whereas the larger To becomes the larger the Mandel parameter becomes, reflecting 
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the large fluctuations implicit in jumps from dark to bright periods. These macroscopic 
fluctuations are manifested in the photon counting distributions studied in detail by Schenzle 
and Brewer (Schenzle and Brewer, 1986) using Bloch equations. They showed that the count 
distribution of photons detected from the strongly allowed transition were Poissonian except 
for an excess of zero counts. In an interval T of order of the lifetime of the shelving state 
one either counts a large number of photons (a bright period) or one counts nothing (a dark 
period). The probability of counting n photons in time T is Poissonian except again for 
an excess of zeros (Schenzle and Brewer, 1986) and is evaluated from the Mandel counting 
formula (or its quantum equivalent derived in (Kelley and Kleiner, 1964)) 

except for n = 0. Here 71 is the decay rate of the strongly-fiuorescing state and r] is the 
detector efficiency. So in an interval T ~ ^2^ where 7^^ is the lifetime of the shelving 
state, we count either n ~ \'Jii]'l2^ ~ 10^ fo^ typical transitions or we find n = 0. Note we 
essentially do not see n = 1, 2, 3, ... as W{n = 1, T = '^2'^) ~ lO^e^^*^ = 0, as shown in Fig. 
6. Schenzle and Brewer (Schenzle and Brewer, 1986) interpret these macroscopic intensity 
fiuctuations in terms of quantum jumps. Imagine the fiuorescent intensity to be jumping 
from dark (off) states to bright (on) states with probability distribution 

PiI) = A6iI) + il-A)5iI-Io) , (41) 

then we find for the probability of counting n photons in a time T 

W{n, T) = A5„,o + (1 - A)^-^^^e-^^<^^ , (42) 

n! 

where H is the counting efficiency. The zero count is W{n = 0,T) = A = 1/3 for saturated 

transitions. The behaviour of Eq. (42) is schematically shown in Fig. 6 where we see the 

excess probability for no counts (dark period) together with a high probability for a large 

number of jumps (bright period). 

In the past two sections we have discussed the initial attempts towards a theoretical 

description of single ion resonance fiuorescence. However, these attempts did not yet yield 



a satisfying approach to the problem as the single system properties described, e.g., by the 
delay function, were deduced from equations of motion describing the whole ensemble. In 
the next section we will now explain and summarize a number of approaches all giving the 
quantum jump approach that allows the most natural description of many properties of 
resonance fluorescence and time evolution of single quantum systems. 

IV. DISCUSSION OF DIFFERENT DERIVATIONS OF THE QUANTUM JUMP 

APPROACH 

A. Quantum jumps 

Prior to the development of quantum jump methods, all investigations of the photon 
statistics started out from the ensemble description via optical Bloch equations, or rate 
equations as presented above, which were used to calculate nonexclusive 'probability densi- 
ties' for the emission of one or several photons at time ti, . . . ,tn in the time interval [0,t]. 
It is important to note that only the probability of emission of any photon was asked for. 
Therefore many more photons might have been emitted in between the times ti. An exam- 
ple for such a function which we discussed in section III is the intensity correlation function 
g'^'^\t) which gives the normalized rate at which one can expect to detect photons (any pho- 
ton rather than the next) at time t when one has been found at t = 0. 
Efforts were made to use nonexclusive "probability densities" to deduce the photon statis- 
tics of the single three-level ion and the aim was to show that a single ion exhibits bright 
and dark periods in its resonance fluorescence on the strong transition (Pegg and Knight, 
1988a; Schenzle and Brewer, 1986). This approach is, however, not very satisfying as it re- 
quires the solution of the full master equation and the inversion of a Laplace transformation. 
Furthermore this approach is very indirect as we first calculate the ensemble properties and 
then try to derive the single particle properties. It would be much more elegant to have a 
method which enables us to calculate the photon statistics of the single ion directly. This 
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was discussed widely at a workshop at NORDITA in Copenhagen in December 1985, follow- 
ing a paper at that meeting by Javanainen (Javanainen, 1986b). This intention was finally 
realized with the development of the quantum jump approach. Its development essentially 
started when Cohen- Tannoudji and Dalibard (1986) and much at the same time ZoUer, 
Marte and Walls (1987) derived the exclusive probability Po{t) that, after an emission at 
time t = 0, no other photon has been emitted in the time interval [0,t] (Cohen- Tannoudji 
and Dalibard, 1986) or the exclusive n-photon probability density P[o,t]it\ti, ■ ■ ■ ,tn) that in 
[0, t] n photons are emitted exactly at the times ti, . . . , t„ (Zoller et ai, 1987) without going 
back to the master equation of the full ensemble. Both quantities are intimately related, as 
the probability density /i(t) for the emission of the first photon after a time t is given by 

h{t) = Poit) , (43) 

and because it turns out that the exclusive n-photon probability density essentially factorizes 
into next photon probability densities /i(t) 

P[0,i](t|tl, . . . ,t„) = Po{t - tn)h{t,, - tn-l) ■ . . . ■ /l(ti) . (44) 

This factorization property was initially assumed and then justified by physical arguments by 
Cohen- Tannoudji and Dalibard in (Cohen- Tannoudji and Dalibard, 1986) while in (Zoller 
et ai, 1987) first the exclusive n-photon probability p[o,t](t|ti, . . . , t„) was calculated and 
subsequently its decomposition into next photon probabilities was derived. Before we discuss 
the approaches of Cohen- Tannoudji and Dalibard (Cohen- Tannoudji and Dalibard, 1986) 
and Zoller, Marte and Walls (Zoller et ai, 1987) we point out that although the exclusive 
(next photon at time t) Ii{t) and nonexclusive distributions (a photon at time t) I{t) are 
very different, they are related by a simple integral equation (Kim et ai, 1987). We have, 
as discussed in Eq. (32) in Section III 

/(t) = /i(t) + /* dt'h{t - t')I{t') , (45) 

Jo 

which becomes especially simple as we saw when one considers the Laplace transform of 

that equation as we have a convolution on the right hand side of Eq. (45). This relationship 
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enables us in principle to obtain the exclusive probability density Ji (t) from the nonexclusive 
quantity I{t). In practice however this is exceedingly difficult to do as one has to know all the 
eigenvalues of the corresponding Bloch equations. Therefore a direct approach is needed. 
The idea put forward in (Zoller et ai, 1987) was to calculate not the complete density 
operator p irrespective of the number of photons that have been emitted but to discriminate 
between density operators corresponding to different numbers of emitted photons in the 
quantized radiation field. The quantity of interest is therefore 

pf)(t)=tr^{IP„p(t)} , (46) 

where p{t) is the density operator of atom and quantized radiation field, trpi-} the partial 
trace over the modes of the quantized radiation field and IP„ the projection operator onto 
the state of the quantized radiation field that contains n photons. This projector is given 
by 

IPn = ;^ E • • • E <A. • • •aLAjO)(0|ak„A„ . . .ak.A. • (47) 

klAl knAn 

This method to calculate the density operator for a given number of photons in the quantized 
radiation field was first used by Mollow (1975) to investigate the resonance fiuorescence 
spectrum of two-level systems. However, as at that time the investigation of single ions was 
completely beyond then-current experimental possibilities, so that he did not draw further 
conclusions from his approach concerning single quantum systems. This was only triggered 
later by the experimental realization of single ions in ion traps. 

In the following we discuss the approach of Zoller et al (1987) for a three-level system 
in V configuration (see Fig. 1) case rather than as in the original paper for the two- level 
case. Following (Mollow, 1975; Blatt et al, 1986) one obtains for Eq. (46) the equations of 
motion 



f/^\t) = -i{H^ffpf{t) - pTit)Hy (48) 

and 
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^A\t) = -^iH^JfP^-\t)-p^-\t)Hy+J:2Tu\0){^\pi-'\t)\^){0\ , (49) 

where the effective Hamilton operator H^ff is given by 

H^ff = -Y.l^HA, + ^ru)\^){^\ + ^m{^\ + \^)mj , (50) 

with the detunings Aj = cjj — un and, Ui the laser frequency, Qi the Rabi frequency and 
2 Tii the Einstein coefficient on the i ^^ 1 transition. It is now important to note that the 
effective Hamilton operator H^jf is a non-Hermitean operator. The real part of —i Heff is 
negative which implies that the trace of the density operator Pj^ {t) decreases in time. This is 
not surprising because p^ (t) describes the conditional time evolution under the assumption 
that no photon has been emitted into the quantized radiation field. The probability that an 
excited atom has not emitted a photon decreases in time and therefore the trace of pV(^) 
describing this probability should decrease in time. This decrease is necessary for the trace 
of the density operator p{t) disregarding the number of emitted photons, 

oo 

T(t,0)p(0) = p(t) = EpJ)(t). (51) 

n=0 

to be preserved under the time evolution. The equations of motion Eqs. (48) and (49) have 
the solution 

pfit) = Sit,to)pfito) , (52) 



where 



and 



with 



S{t,to)pf{to) = e-^^^//(*-*")p5')(to)e^^«^//(*-*") (53) 



p^^\t)= fdt'Sit,t')Rpt'\t') (54) 

Jo 



i?p!rHt) = ^2rdo)(^ip!rHt)K)(oi. (55) 



j=l 
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From this result it is now possible to deduce the probability that exactly n photons have 
been emitted in the time interval [0,t]. The probability that no photon has been found 
should then be given by 

Po(t)=tr^{5(t,0)pA(0)} , (56) 

and 

Pn{t) = f dtn... r dtitrA{S{t,tn)R...RS{h,{))pAm (57) 

JO Jo 

is the probability that exactly n photons have been emitted. Zoller, Marte and Walls then 
realized that the structure of these expressions coincides with that derived from an abstract 
theory of continuous measurement constructed by Srinivas and Davies (Davies, 1969; Davies, 
1970; Davies, 1971; Davies, 1976; Srinivas and Davies, 1981; Srinivas and Davies, 1982). This 
theory supports the interpretation of 

p„(ti, . . . , t„) = trA{S{t, tn)R . . . RS{h, O)pa(O)} (58) 

as the probability density that exactly n photons have been emitted at times ti, . . . , t„ and 
no photons in between. From the general theory of measurement, they interpreted the 
quantity 

p ., , . _ tr^{i?g(ti,to)i?T(to,0)p^(0)} 

"^'^'^ " '°^ - tr^{i?T(to,0)p^(0)} ('^^ 

as the probability density that after an emission of a photon at time to the next photon 
will be emitted at ti. It should be stressed that although in (Zoller et ai, 1987) the super- 
operator T{t, 0) is identified with the time evolution of the ensemble irrespective of how 
many photons have been emitted in [0, t] for the following it should always be chosen to be 
the time evolution if a given number of emissions have taken place at the times ti, . . . ,t„. 
Assuming this (as is also implicitly done later in (Zoller et ai, 1987)) one obtains for the 
probability density that photons are emitted exactly at times ti, . . . , t„ 
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P[o,t]iti, ...,tn) = trA{S{t, tn)R . . . RS{tu O)pa(O)} 

trA{Sit,tn)R...RSih,0)pAm trA{RS{t,,0)pAm 
trA{RS{t, tn-i) . . . RS{tu O)pa(O)} ■ ■ ■ tTAipAiO)} 

= Po(t,t„)/i(t„,t„_i).../i(ti,0) . (60) 

Here we have factorized p[o,t](ti, . . . , t„) into products of Ii{ti, ti^i). In principle these func- 
tions can depend on the atomic state at time t;_i (after the emission). However, in most cases 
this state will be the ground state of the system and will be the same after each emission. 
Having found that the knowledge of Po{t) is sufficient (/i(t) can be obtained via Eq. (43)), 
Zoller et al then continue to discuss the photon statistics of the three-level V system. The 
results found in (Zoller et al, 1987) may also be used to implement a simulation approach 
for the time evolution of a single three-level system (Dalibard et al., 1992; Dum, Zoller and 
Ritsch, 1992; Hegerfeldt and Wilser, 1991). However, the application of the quantum jump 
approach in numerical simulations will be discussed later in this section. 

The approach of Zoller, Marte and Walls already reveals many features of the quantum 
jump approach. However there is a slight complication in their approach as they rely on the 
abstract theory of continuous measurement of Srinivas and Davies to give interpretations 
to their expressions Eqs. (58)-(59). The reason that they need the support of the theory 
of Srinivas and Davies is that they never talk about the way the photons are measured. In 
fact only the emission of photons is mentioned and not the detection of photons. From a 
quantum mechanical point of view, however, one has to be very careful, as the emission of a 
photon is not well defined. It is the detection of a photon in the radiation field which is a real 
event. Of course the treatment in (Zoller et al., 1987) already implies some properties of the 
measurement process, e.g., they implicitly assume time resolved photon counting. However, 
no explicit treatment of such measurements was given in (Zoller et al., 1987). This problem 
was then addressed in greater detail by several authors and in the following we discuss these 
ideas. 

The first approach to include the result of quantum mechanical measurements into their 
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calculation explicitly was given by Porrati and Putterman (1989). They (Porrati and Put- 
terman, 1987) as well as others (Pegg and Knight, 1988a) noticed that the failure to detect 
a photon in a measurement leads to a state reduction, as information is gained through 
that null measurement. Essentially we can be increasingly confident that the ion is in a 
non radiating state (examples of this will be shown in Section V). Porrati and Putterman 
assume that at some large time t a measurement on the whole quantized radiation field is 
performed. Assuming the result of this measurement is the detection of no photons, they 
calculate all Heisenberg operators at that time projected onto the null-photon subspace of 
the complete Hilbert space, i.e., operators of the form 

(IPo6AlPo)(t) . (61) 

Although not mentioned explicitly in (Porrati and Putterman, 1987; Porrati and Putter- 
man, 1989), the calculation of this operator turns out to be closely related to the projector 
formalism (Agarwal, 1974; Haake, 1973; Nakajima, 1958; Zwanzig, 1960), a connection that 
was elaborated on by Reibold (Reibold, 1992). Although their approach can in principle lead 
to the quantum jump method, there are some conceptual problems in the actual execution 
of the use of the null- measurements idea. The main problem is that Porrati and Putter- 
man only talk about a single measurement at a large time t performed on the complete 
quantized radiation field. This does not seem to be a very realistic model of measurements 
performed by a broadband counter, which informs us immediately whether he has detected 
a photon or not. Also the calculation of the state after the detection of a photon was not 
elaborated on in (Porrati and Putterman, 1987; Porrati and Putterman, 1989), where it was 
merely stated that the system is reset back to its ground state on photodetection which 
is of course a physically correct picture. These conceptual concerns to this approach were 
later addressed in the work of Hegerfeldt and Wilser (Hegerfeldt and Wilser, 1991; Wilser, 
1991), of Carmichael (Carmichael, 1993a) and of Dalibard, Castin and M0lmer(Dalibard et 
ai, 1992; M0lmer et ai, 1993; M0lmer, 1994) and in the following we give a more detailed 
account of their approach. 
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We will follow closely the presentation given by the Hegerfeldt group (Hegerfeldt and 
Wilser, 1991; Wilser, 1991) as it directly leads to the delay function that was also found 
in earlier papers ( Cohen- Tannoudji and Dalibard, 1986; Zoller et ai, 1987). The physical 
ideas, however, are very similar to those presented elsewhere (Dalibard et ai, 1992; M0lmer 
et ai, 1993; M0lmer, 1994). We treat the same three level system as in the discussion of 
Zoller et al (1987). 

In Hegerfeldt and Wilser (1991) and Wilser (1991) (as well as in (Dalibard et al, 1992; 
M0lmer et al., 1993; M0lmer, 1994)), the following simple model of how the photons are 
detected was proposed. It was assumed that the radiating ion is surrounded by a 47r- 
photodetector that detects photons irrespectively of their frequency and that the efficiency 
of the detector is unity. Efficiencies less than unity may be treated in a mathematically 
slightly different way using the same physical ideas (Plenio, 1994; Hegerfeldt and Plenio, 
1996) and leads to a natural connection between the next photon probability density /i(t) 
and the any photon rate (intensity correlation function) g'^'^\t) (Plenio, 1994; Kim et al., 
1987). We will return to this point later. As truly continuous measurements in quantum 
mechanics are not possible without freezing the time evolution of the system through the 
Zeno effect (Mahler and WeberruB, 1995; Misra and Sudarshan, 1977; Reibold, 1993), it 
was instead assumed that measurements are performed in rapid succession where the time 
difference At between successive measurements should be much larger than the correlation 
time of the quantized radiation field. This means that 

At > uo^^ . (62) 

If successive measurements are more frequent than ojiq we enter the regime of the quantum 
Zeno effect and we significantly inhibit the possibility of spontaneous emissions (Reibold, 
1993). On the other hand At should be much smaller than all time constants of the atomic 
time evolution to ensure that one finds the photons one by one and because we want deter- 
mine the time evolution using perturbation theory. Therefore 

r,::\Ar\f]ri«At . (63) 
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For optical transitions it is easy to satisfy both inequalities Eqs. (62) and (63) simultane- 
ously. 

Now the density operator at time t under the condition that no photons have been 
detected in all measurements which took place at time si,...,s„ has to be calculated. 
Although the result of each measurement was negative, in the sense that no photon was 
found, this still has an impact on the wavefunction of the system, as it represents an increase 
of knowledge about the system (Dicke, 1981; Porrati and Putterman, 1987; Pegg and Knight, 
1988a). Using the projection operator IPq onto the vacuum state of the quantized radiation 
field and the time evolution operator U{t, t') of system and radiation field in a suitable 
interaction picture, we find 

\i,{sn)) = IPof/(s„,s„_i)IPo... IPof/(si,0)|z^(0)) . (64) 

as, after each measurement which has failed to detect a photon, we have to project the quan- 
tized radiation field onto the vacuum state according to the von Neumann-Liiders postulate 
(Liiders, 1951; von Neumann, 1955). As Sj — Sj_i obeys Eqs. (62) and (63), it is possible 
to calculate the time evolution operator f/(sj, Sj-i) in second order perturbation theory to 
obtain 

IPo?7(s„s,_i)IPo ^ {l-'-H,ff{si.i){s, - s._i))IPo (65) 

with the effective Hamilton operator Eq. (50). In the quantum jump method presented in 
(Dalibard et ai, 1992; M0lmer et ai, 1993; M0lmer, 1994) the Weisskopf-Wigner approach 
(Weisskopf and Wigner, 1930) was used to find a formula equivalent to Eq. (65). Inserting 
into Eq. (64) and going over from a coarse grained time scale to a continuous time, we 
obtain for the atomic part of the wavefunction (the radiation field is in its vacuum state) 
where no photons have been detected in all measurements in the interval [0,t] 

|^f(t)) = e-^^//*/^|^^(0)). (66) 

One should note that the effective time evolution does not preserve the norm of the state 
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and that it maps pure states onto pure states. In fact the square of the norm of Eq. (66) is 
just the delay function 

Po{t) = {^f{tMf{t)) (67) 

and coincides with Po{t) given in Eq. (56). The delay function Po{t) will become important 
in applications of the method in simulations (Dalibard et ai, 1992; Dum, Zoller and Ritsch, 
1992). It should be noted here in passing that if we consider the normalized version of the 
time evolution Eq. (66) for a two-level system, then one finds that it is identical to the 
time evolution according to the neoclassical radiation theory of Jaynes (Bouwmeester et ai, 
1994). The reason for this is claimed by Bouwmeester et al to rest on the fact that H^/f 
contains all contributions from virtual photons (i.e. all radiation reaction terms) but does 
not include the real photons as their detection leads to state reduction according to the 
projection postulate. However, neoclassical theory predicts quantum beats from a three- 
level system in A configuration while it is easily seen that an analysis of the problem using 
the quantum jump approach does not predict quantum beats; a result in accordance with 
experiment Milonni (1976) and references therein. 

Eventually the photo detector will find a photon and the state after this detection can 
be determined by the projection postulate. We write the state after the detection of a 
photon as a density operator, as the state after the emission can be a mixture (e.g. as in 
the three-level A configuration (Javanainen, 1992; Hegerfeldt, 1993; Hegerfeldt and Plenio, 
1995a; Hegerfeldt and Sondermann, 1996)) although in the case of the three-level V system 
it is not: 

PR{Sn) = (i- IPo)f/(Sn,Sn-l)IPoP?^(s„_l)IPof/(Sn,S„_l)(i- IPq) . (68) 

At this point the additional assumption is made that the photo detector absorbs the photon 
(which it does in reality) and that the state after the detection is simply obtained by remov- 
ing the photons from the radiation field. This can be done by tracing over the quantized 
radiation field and multiplying with IPq 
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PR = trpipRiSn)} ® IPo 
2 
= J22T,,\0){l\pfiSn-l)\t)mSn-Sn-l) • (69) 



It should be noted that the assumption that the state after the detection of a photon in 
the counter is given by Eq. (69) is not included in the projection postulate but enters as 
a physically justified additional assumption. However using slightly different mathematical 
methods it is possible to show that the procedure Eq. (69) is not really necessary. Ideal 
quantum mechanical measurements where the photon is not absorbed lead to the same 
results (Plenio, 1994; Hegerfeldt and Plenio, 1996). This is a consequence of the intuitively 
obvious fact that in free space photons emitted by the system will never return to it and is 
implicit in the treatment of Zoller et al (1987). 

A different approach towards the quantum jump method was presented earlier by Carmichael 
and coworkers (Carmichael et al., 1989; Carmichael, 1993a). They derive the quantum 
jump method from a discussion of photoelectron counting distributions that are found in 
experiments. A quantum mechanical theory for photoelectron counting distributions was 
developed in 1964 by Kelley and Kleiner (Kelley and Kleiner, 1964) who derived the quantum 
mechanical expressions for nonexclusive multicoincidence rates. For the probability to have 
n photoelectron counts in the time interval [t, t + T] , they find 

where ^ is the product of detector efficiency and a factor to convert field intensity into 
photon flux. The notation (: ... :) means that all operators have to be normally ordered 
and time ordered in such a way that times decrease from the centre towards the left and 
right. Expanding Eq. (70) one can write p{n,t,T) as a complicated series of integrals over 
the nonexclusive multicoincidence rates 

/(ti, . . . ,U = r{E^'\h) . . . M-\tjM+\t„,) . . . M^\t,)) , (71) 

which gives the rate for the joint detection of photons at times ti, . . . , tm- It is a nonexclusive 
rate as there may be more detections in between the times ti, . . . ,tm- That these possible 
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events are included in Eq. (71) is obvious, as the Heisenberg operators are calculated with 
respect to the total Hamiltonian of the system which describes a time evolution in which 
arbitrarily many photons may be created. The analysis of the photon statistics by means 
of the coincidence rate Eq. (71) has long been the standard way of investigation. It was, 
however, realized that this is not the only possibility, and for certain problems it is not even 
the most natural way. In fact Eq. (70) , for example, can be expressed very easily by the 
exclusive probability density to find photon counts at exactly the times ti, . . . , t„ and at no 
other time in [t,t + T]. One finds for this the expression 

dtn J dtn-l ... j dtiP[t,t+T]{tl, ...,tn) • (72) 

Carmichael and coworkers (Carmichael et ai, 1989; Carmichael, 1993a) then undertook the 
step to express the exclusive probability density in terms of the intensity operators. They 
find 

P[t,t+T]{ti,...,tn) = Y.^^ <•••/ dt[{:I{t',)...I{t[)I{tJ...I{t^):) , (73) 

where 

i{t)=^E^~\t)E'^+\t). (74) 

This can be checked by inserting Eqs. (73) and (74) into Eq. (72) and showing that the 
result coincides with Eq. (70) (Saleh, 1978; Stratonovitch, 1963). The aim now is to rewrite 
Eq. (73) in terms of (super)-operators that only act in the atomic space, as these are much 
easier to handle than the Heisenberg-operators I(t). It turns out that the resulting equations 
are quite simple. To this end it is important to note that the electric field operator in the 
Heisenberg picture can be decomposed into a free-field part and a source field part 

E^-\t) = Ey\u) + Ei-\U) . (75) 

In the Markov approximation the free field commutes with all electric field operators at 
earlier times Cohen-Tannoudji, J. Dupont-Roc, and G. Grynberg (1992) and when acting 
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onto the vacuum state it vanishes. It is therefore possible to replace the intensity operators 
/(tj) for i = 1, . . . ,m by E^~\ti)E^g'^\t,i) in Eq. (73). All other intensity operators remain 
unchanged. Carrying out the time ordering in Eq. (73) explicitly and after tedious calcu- 
lations (for details we refer the reader to (Carmichael et ai, 1989; Carmichael, 1993a)) one 
obtains for an overall counter efficiency ^ 

P[M+T](ti, . . . , t„) = ^trA{e(^^^^)(*+^-*-)7^ . . . 7^e(^-«^)(*-*)/.(t)} , (76) 

where C is the superoperator given by 

Cp = -^iH,jfiOp-pHeffm , (77) 

where Hg.ff{C,) is obtained from H^ff for perfect efficiency ,^ = 1 by substituting F -^ (1— ,^)r. 
TZ is the reset operator giving the state after a photon has been detected. Assuming unit 
efficiency (.^ = 1) of the detection process we recover Eq. (50). 

So far we have discussed a number of approaches to the quantum jump description of dissipa- 
tion. These approaches can be formulated somewhat differently in the language of quantum 
stochastic differential equations (Gardiner, 1992). This formulation is certainly rather formal 
at first glance, but it has the advantage that certain operations where one uses the Markov 
approximation become simpler. On the other hand one has to use the somewhat unintuitive 
Ito formalism (Gardiner, 1992; Gardiner et al, 1992) and a more physically oriented deriva- 
tion would sometimes be helpful for the interpretation of the occurring equations. 
To illustrate the idea of this formalism we consider a laser-driven two-level atom in a quan- 
tized radiation field which is in the vacuum state. We follow the description in Sondermann 
(1995b). The Hamilton operator in a suitable interaction picture is given by 

H = -^Ai|l)(l| + ^(|0)(1| + |1)(0|) + E(^%kA^ioakAe-*^"''^-'^i«)* + h.c.) 

^ kA 

= Ha + aioDioEW(t) + DloE(-)(t)aoi , (78) 

where aij = \i){j\ is an operator annihilating an electron in level j and creating an electron 
in level i and where 
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/ 2 \ \ (1/2) 

DioEW(t) = Dio E^^ (12^) ^kA^ioakAe-^(-----)* (79) 

is the interaction energy between the electric-field operator -E^(t) in the Schrodinger picture 
(or more precisely in the chosen interaction picture) and the atomic dipole moment D21 of 
the transition. The time discretized Schrodinger equation then reads 

ihA\ilj{t)) = H\ilj{t))At 

= {HaM + aioAAt(t) - a^^AA{t)}\i,{t)) , (80) 



AA{t) = j^ dA{t) = j^ dt Dio E ^^ I y^ ) ekA^ioakAe-^-- --"')* 



where 

(1/2) 

kA V2eo/^\/J 

We assume in the following that At ^ ^h)^, which is crucial for us to be able to perform 
the Markov approximation. The idea is now to perform the Markovian limit directly in 
the Schrodinger equation instead of performing this limit on the results. This is the step 
where we have to introduce the notion of quantum stochastic differential equations, as in 
performing this limit we cannot subsequently interpret the resultant Schrodinger equation 
as an ordinary differential equation anymore (Gardiner, 1992; Gardiner et ai, 1992). Under 
the Markov assumption we have 

if It - si > At 

[AA{t),AA^{s)] = { '- , (82) 

I rAt else 
which in the limit At ^ results in 

1 ift^s 

[dA{t),dA^is)] = I . (83) 

I Flit else 

We also need to know that if we assume that the initial state of the quantized radiation field 
is the vacuum, then 



''1 



dA{t)dA{t) = = dA\t)dA^{t) , (84) 
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where this equation is defined in the mean-square topology sense, i.e., in brief one apphes 
both sides on an initial vector and takes the absolute square of the result afterwards (Gar- 
diner, 1992; Gardiner et al, 1992). Taking the limit At ^ in Eq. (81) we have assumed 
the ordinary rules of calculus, and therefore generated a stochastic differential equation in 
the sense of Stratonovitch 

iM\ij{t))\s = {HAdt + aiodA^t) - aQidA{t)}\ilj{t)) . (85) 

As a Stratonovitch equation is not easy to integrate, we would like to transform it to an 
Ito-form using the rules Eqs. (82)-(84) (Gardiner, 1992; Gardiner et al, 1992). We then 
find 

iMlipit))]! = {HAdt + dA\t)aoi - Tandt}\^lj{t)) , (86) 

where the Taudt arises from a dA{t)dA'^{t) contribution. As the dA{t) commute with all 
earlier dA{s) upon which {ipit)) depends, it can be commuted to the right until it operates on 
the initial state and therefore the vacuum. Therefore the contribution of the dA{t) vanishes 
and only the dA'^it) contribution survives. 

What we are in fact interested in here is the rederivation of the quantum jump approach. 
Therefore we are interested in the time evolution when no photon is present in the field, i.e., 
we are interested in the state vector 

IPol^(t)) = \m)o , (87) 

where IPq is the projector onto the vacuum state of the quantized radiation field. We find 

d\m)o = Pod\m) 

= i-tHA/h - Tan)Pom))dt • (88) 

The (iA^(t) now vanishes because acting on a vacuum state to its left it gives zero contri- 
bution. The norm of the conditional state vector |-?/'(t))o is just the probability to find no 
photon until t if there was no photon at t = 0. This is just the reduced time evolution found 
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in the previously discussed derivations other approaches, too. The probabihty density Ii{t) 
for an emission at time t is just the rate of decrease of the norm of the emission free time 
evolution, i.e., 

do{m\m)o\i = (doimi) \m)o +o{m\ id\m)o) + (doimi) id\m)o) 

= o{^it)\iHA/h - ran|^(t))o + o(^(t)| - iHa/Ti - raii|^(t))o 

= -2ro(^(t)kii|^(t))o . (89) 

Therefore using this formulation we have recovered the quantum jump approach; we observe 
that this formalism, although not delivering new insights into physics different from those 
from previous derivations of the quantum jump approach, is very elegant from a formal 
point of view. To understand the formalism a little better, we now show how one may 
obtain Eq. (88) without referring to the formalism of stochastic differential equations (ZoUer 
and Gardiner, 1995). We consider a finite time step for the state vector IPo|'?/'(t)) using the 
Hamilton operator Eq. (78) and first order perturbation theory. We obtain 

AlPolV^(t)) = \P^{-iHAin/\t - aioAA(t))|^(t)) 

= IPof/(t, 0)f/(^)(t, 0){-iHA/hAt - aioAA{t))U{t, 0)U'^^\t, 0)|^(t)) 

= -tHA/nAtPoim) + IPof/(t,0)aio {AA{t) - Taoi{t)) |^(0)) , (90) 

where in the last line we have used the well known expression for the Heisenberg operator 
of the electric field operator which can be written as the free field contribution and a source 
term (the dipole of the atom radiates the outgoing field) (Loudon, 1983). Note that o-oi(^) 
is now a Heisenberg operator. Eliminating A{t) in the last row of Eq. (90), as it operates 
on the initial vacuum state, we obtain 



^5 



AlPol^(t)) = At i-tHA/n - Tan) Poim) ■ (91) 

Now we may easily perform the limit At ^ to obtain the same result as in Eq. (88), 
however, without the explicit use of the quantum stochastic differential calculus. 
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B. Quantum state diffusion and other approaches to single system dynamics 

So far we have discussed the quantum jump approach for the description of single ra- 
diating quantum systems. The main ingredient in the derivation was the assumption of 
time resolved photon counting measurements on the quantized radiation field. The re- 
sulting time evolution could be divided into a coherent time evolution governed by a non 
Hermitean Hamilton operator which is interrupted by instantaneous jumps caused by the 
detection of a photon and the consequent gain in knowledge about the system. One could 
ask whether this description is unique, that is, it represents the only possibility. From the 
emphasis we put on the importance of the measurement process in the derivation of the 
quantum jump approach one can already guess that other measurement prescriptions will 
yield different kind of quantum trajectories. In the following we will discuss an important 
example, quantum state diffusion, (Gisin and Percival, 1992a) of a different kind of quan- 
tum trajectories which in fact can be derived from a very important measurement method 
in Quantum Optics, namely, the balanced heterodyne detection. Before we show the con- 
nection of quantum state diffusion to balanced heterodyne detection let us point out that 
quantum state diffusion was originally derived independently from a measurement context. 
Steps in this direction were made when several authors became interested in alternative ver- 
sions of quantum mechanics (Pearle, 1976; Ghirardi et ai, 1986; Ghirardi et al, 1990; Diosi, 
1988; Diosi, 1989) and the investigation of the wavefunction function collapse, i.e., the pro- 
jection postulate (Gisin, 1984; Gisin, 1989). In these investigations stochastic differential 
equations for the time evolution of the state vector of the system were studied. Again there 
is a multitude of possible equations; however, Gisin and Percival (1992a) provided a natural 
symmetry condition under which it is possible to derive a unique diffusion equation which is 
referred to as the quantum state diffusion model (QSD). Given a Bloch equation in Linblad 
form (Lindblad, 1976) 

P = - ^ [Hsys, p]+J2(^ LrnpL^m - ^m^mP " pL'^mLm) , (92) 
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with the system Hamiltonian Hgys-, and the Lindblad operators L^, the quantum state 
diffusion equation for the state vector is 

+ J2{Lm-{Lm))\^)dU ■ (93) 

m 

The d^m represent independent complex normahzed Wiener processes whose averages, de- 
noted by M(. . . ), satisfy 

M{dU) = , 
M{Re{dU)Re{dU) = M{Im{dU)Im{din)) = Smndt , 

M{Re{dU)Im{d^n)) = . (94) 

Equation (93) has to be interpreted as an Ito stochastic differential equation, see , for 
example, (Stratonovitch, 1963; Gardiner, 1992). It is easy to check that averaging Eq. (93) 
over the stochastic Wiener process yields the density operator equation Eq. (92) and that 
therefore (in the mean) normalization is preserved. For numerical studies often a somewhat 
simpler equation is used that does not preserve normalization even under the mean. This is 
given by 

|#) = - ^H\^)dt + ^ {2{Ll)Lm - LIlJi \^)dt + Y, L„,\i,)dU . (95) 

''' m m 

It should be noted that Eq. (94) is a nonlinear equation as it also depends on the expectation 
values of the Lindblad operators L^- This makes the analytical treatment of this equation 
very difficult and there are only a few cases for which an analytical solution is known (Gisin, 
1984; Gisin, 1989; Salama and Gisin, 1993; Wiseman and Milburn, 1994; Carmichael, 1994). 
However it was found by Goetsch, Graham and Haake (Goetsch and Graham, 1993; Goetsch 
and Graham, 1994; Goetsch et ai, 1995) that it is possible to find linear stochastic differen- 
tial equations which also reproduce the ensemble average. Stochastic differential equations 
for the wavefunction have also been derived by Barchielli (Barchielli, 1986; Barchielli and 
Belavkin, 1991; Barchielli, 1993) (see also (Zoller and Gardiner, 1995) for good summary 

36 



of these approaches) from a more abstract mathematical point of view. The approach of 
Barchielh also gives a common mathematical basis for both diffusion and jump processes. 

However, we do not intend to elaborate further on the mathematical side of the theory. 
Instead we would like to show that it is possible to derive QSD from the quantum jump 
approach in a certain limiting case, i.e., the case of infinitely many jumps where each jump 
has an infinitesimal impact on the wavefunction. In fact it turns out that QSD can be related 
to an explicit and well known physical measurement process in quantum optics, namely, the 
method of balanced heterodyne detection, see , for example, (Castin et ai, 1992; Wiseman 
and Milburn, 1993a; Wiseman and Milburn, 1993b; Wiseman, 1996; Carmichael, 1993a; 
M0lmer, 1994; Knight and Garraway, 1996). In the following we would like to show this 
explicitly for the specific example of a decaying cavity and we follow a similar path to that 
used in the approach of Garraway and Knight (Gastin et ai, 1992; Garraway and Knight, 
1994b; Knight and Garraway, 1996). To be specific we will illustrate the method for the 
case of balanced heterodyne detection of the output of an undriven optical cavity. We have 
in mind the situation given in Fig. 7. 

The left hand cavity A (with mode operators acav) is the source a weak output field 
(mode operators CkA) which we want to analyze, while the lower cavity (mode operators 
bloc) is assumed to be in a coherent state \f3) with a very large amplitude f3 for all times so 
that the radiated field (with mode operators 6kA) of that cavity is very large. The Hamilton 
operator describing this situation is 

H = huJcavCtlav'^cav + Yl ^^kAOk^akA + Yl {^^akXalavO-kX + ^-cj 

kA kA 

+ huJioch\ochloc + Y ^^kA^A^kA + Y {^^/kA&L^kA + h.C.^ . (96) 

kA kA 

where the (7kA and /kA are the coupling constants between the cavity and the outside world 
and Ucav and uioc the frequencies of the cavity A and the local oscillator cavity respec- 
tively. The action of the beamsplitter is to mix the two incoming modes. Assuming a 50% 
beamsplitter we find for the new mode operators (Loudon and Knight, 1987) 
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C?kA = -/^(-OkA + &kA) • (97) 

Now going over to an interaction picture with respect to 

Ho = hUcavalavO'cav + Y^ ^t^kAOk^akA + ^^loMoMoc + Yl ^^kA&L^kA (98) 

kA kA 

and subsequently changing the basis via a displacement operator such that the initial state of 
the local oscillator is the vacuum (Mollow, 1975; Pegg, 1980) we obtain using Q = uioc — oJcav 
ih 



^ = E ^ (^kAaLe-^^* + /kA(&L + Pi) CkAe 

kA V^ 

ih 



+ E ^ (-^kAaLe-^™ + /kA(&L + n) 4Ae"^(---'-)* + h.c. . (99) 

kA v2 ^ 

Now applying the methods that we used to derive the quantum jump approach, we easily 
obtain the two jump operators 



^d = ^ (-VT^Oca^e*™ + ^J^cP) , (100) 

where "-fcav and 7zoc are the decay rates of the cavity and the local oscillator. Within a short 
time interval At, i.e., such that {oJcav — ^ioc)^t ^ Ij we will count on average 



i-PM - % (i + /If (-«.)) (101) 



counts in mode c and 



(^j^.>-%(i-/;§(.o,)) . (102) 

counts in mode d where 

(^0):=(«Le^'^ + «ca.e-*^). (103) 

These are average values around which the actual number of counts in the two counters 
fluctuates. We can approximate this number of counts m{t) by the stochastic process 
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[t) = (jV)At + {J^jy/^AW (104) 



such that {{AWy) = At. For the powers of the jump operators we find 






As we normahze after each emission the prefactors are not really important and we can 
divide the jump operators by these. One should note that the phase of these prefactors is 
fixed due to the fact that in the limit of infinite (3 the jump operators have to become the 
unit operator. Therefore there is no freedom in the choice of the sign of the prefactors. It is 
now easy to derive the effective Hamilton operator -f^e// for which we find 

H,sf = -th^ia^a + ^^P^) . (106) 

Using this together with Eq. (105) we obtain 



m + At)) = (1 - ^iJe//At + e-*™a,,,(27ea„(xnt) + ^^(Al^i - AW2))m)) ■ (107) 

Adding together the two Wiener noises AW = (AWi — AW2) / V^, taking the limits At -^ dt 
and AW — > dW and defining 

di = e-'^'dW , (108) 

we obtain after dropping a counter-rotating term of the form e~^*^* 

l#) = 



^ Icav 



(109) 



This is the unnormalized diffusion equation given , for example, by Gisin and Percival (Gisin 
and Percival, 1992a). One should note that if we had considered homodyne detection, i.e., 
the case fi = 0, then we would have found a different diffusion equation, as there would not 
have been a counter-rotating term which we could have dropped. Therefore an additional 
term in Eq. (109) would appear (Carmichael, 1993a; M0lmer, 1994). 
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To yield the normalized equations for QSD as they are given by Gisin and Percival 
(1992a) we have to normalize the wavefunction and we have to include a stochastic phase 
factor a{t) into the wavefunction (Garraway and Knight, 1994b) , i.e., we look for a diffusion 
equation for 

\m) = ^7^ ■ (110) 

The reason we have to include this seemingly unmotivated phase factor is that in the deriva- 
tion of the QSD equation in (Gisin and Percival, 1992a), a term is added to the diffusion 
equation to give it the simplest possible form. This term in fact gives rise to a random phase 
change. To yield QSD we choose a(t) as 

"W = Y(aca^)c?^- Y(aL)c?r • (111) 

This choice has the effect of removing dC,* that would appear in the diffusion equation of the 
normalized wavefunction without the additional phase factor. Using Eq. (Ill) in Eq. (109) 
and assuming dC,dC,* = dt we finally obtain 

7 11 

l#) = -2^Lv(^cav + l{acav{alav) + ^^^l^vf^cav) " ^{o-lav) {acav))\i') dt 

+ y/l{acav - {acav))di . (112) 

We have therefore shown that the quantum state diffusion equation Eq. (112) (or Eq. 
(93) ) can be regarded as a limiting case of the quantum jump approach. In Section V 
we will illustrate the transition from the quantum jump behaviour to the quantum state 
diffusion behaviour which takes place when we increase the amplitude of the local oscillator 
(Granzow, 1996). It should be noted that it is also possible to obtain a jump-like behaviour 
from quantum state diffusion equations. However, this procedure is much less satisfying than 
the above derivation of quantum state diffusion from quantum jumps. The reason is that 
one has to modify the quantum state diffusion equation by adding an additional operator, 
the localisation operator. The amplitude with which this localisation operator appears in 
the equations is arbitrary and has to be adjusted according to the experimental situation 
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(Gisin et ai, 1993). This is not particularly satisfying at least in cases in which we deal with 
single ion resonance fluorescence. Here the quantum jump approach appears to be much 
more natural. In fact it can be shown that one can not associate the jumps occurring in the 
quantum state diffusion picture with photon emissions as such an interpretation can lead to 
more than one emission from an undriven two- level system (Granzow, 1996). Taking these 
considerations into account one could be tempted to say that the quantum jump approach is 
more fundamental then quantum state diffusion. However, both approaches have the same 
justification as they were both derived from a particular measurement situation. Depending 
on the experimental situation and the measurement scheme employed we have to choose 
either the quantum jump approach or the quantum state diffusion model to obtain the correct 
description of the experimental situation. Quantum state diffusion was, as mentioned before, 
not originally introduced to describe a specific experimental situation. It was rather seen as 
an attempt to formulate alternative versions of quantum mechanics and there are attempts to 
derive diffusion equations from fundamental ideas such as , for example, decoherence induced 
by gravitational fluctuations (Percival, 1994b; Percival and Strunz, 1996; Percival, 1995). 
Although the quantum state diffusion model can not be regarded as the proper description 
of quantum jumps in single photon counting experiments but rather as the description of 
heterodyne detection it is nevertheless useful in the investigation of single system behaviour. 
Interesting phenomena such as localisation in phase and position space (Gisin and Percival, 
1993a; Gisin and Percival, 1993b; Herkommer et al, 1996; Percival, 1994a) are found. These 
can be used to improve the performance of simulation procedures using a "moving basis" 
approach (Schack et ai, 1995; Schack et ai, 1996a) where only time dependent subset of all 
basis states is used in the simulation. A similar method is also possible for a variant of the 
quantum jump approach (Holland et al, 1996). 

So far, we have discussed the evolution of open systems, that is of microsystems in 
contact with Markovian reservoirs such as the bath of vacuum field modes responsible for 
spontaneous emission. The quantum jump concept within an open system context has to 
do with the gain in information about the microsystem which is accessible from the record 
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available in the dissipative environment. Such jump processes do not require an extension 
or modification of conventional quantum mechanics, and we refer to these as "extrinsic" 
jumps. A very different jump mechanism has been studied by a number of authors (Diosi, 
1989; Ghirardi et al, 1986; Ghirardi et ai, 1990; Percival, 1994b; Percival and Strunz, 
1996). In these approaches the Schrodinger equation is modified in such a way that quantum 
coherences are automatically destroyed in a closed system by an intrinsic stochastic jump 
mechanism. This should be distinguished from the extrinsic mechanisms we are concerned 
with in the bulk of this review. 

To see how an intrinsic jump mechanism works, we need a concrete realisation which 
we can apply to a specific time evolution. Milburn has proposed just such a realisation 
(Milburn, 1991), in which standard quantum mechanics is modified in a simple way to 
generate intrinsic decoherence. He assumes that on sufficiently short time steps, the system 
does not evolve continuously under normal unitary evolution, but rather in a stochastic 
sequence of identical unitary transformations. This assumption leads to a modification 
of the Schrodinger equation which contains a term responsible for the decay of quantum 
coherence in the energy eigenstate basis, without the intervention of a reservoir and therefore 
without the usual energy dissipation associated with normal decay (Moya-Cessa et ai, 1993). 
The decay is entirely of phase-dependence only, akin to the dephasing decay of coherences 
produced by impact-theory collisions or by fiuctuations in the phase of a laser in laser 
spectroscopy, but here of intrinsic origin. 

It is interesting to apply Milburn's model of intrinsic decoherence to a problem of dynam- 
ical evolution: that is, the interaction of two subsystems and the coherences which establish 
themselves as a consequence of their interaction. In (Moya-Cessa et ai, 1993) the interac- 
tion between a single two-level atom and a quantized cavity mode was considered and shown 
how the intrinsic decoherence affects the long-time coherence characteristics of the entangled 
atom-field system. In particular it could be shown how the revivals (a signature of long- 
time coherence) are removed by this intrinsic decoherence. The quantum jump approach as 
we have discussed it so far only treats systems (atoms) interacting with a Markovian bath 
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(the quantized multimode radiation field). However, one might be interested to apply the 
quantum jump approach to non Markovian interactions. Examples are electrons interacting 
with phonons or, in quantum optics, an atom in a cavity interacting with a mode which loses 
photons to the outside world (Garraway, Knight and Steinbach, 1995). The second example 
already suggests a possible way one could model such systems. Here the atom sees a cavity 
mode with finite width, i.e., a spectral function which is not fiat but a Lorentzian (Piraux 
et al, 1990 and references therein. However, one does not need to solve a non Markovian 
master equation, as the width of the mode is produced by its coupling to the outside world. 
Taking this coupling explicitly into account, by describing a coupled atom-cavity field mode 
with a dissipative field coupling to the environment, one again obtains a Markovian master 
equation. This is also the recipe for the treatment of an interaction with a bath with a 
general spectral function R{u) (Imamoglu, 1994). One has to decompose R{uj) into a sum 
(or integral) of Lorentzians with positive weights. Each Lorentzian can then be modelled 
by a mode interacting with both the system and a Markovian reservoir. This method is 
practical only if the number of additional modes that one has to take into account is not 
too large. One should also note that in this case the meaning of a jump in the simulation 
can become obscure, as the excitation of the system is transferred to the Markovian bath in 
two steps via the additional mode (Garraway, Knight and Steinbach, 1995). However, if one 
is only interested in a simulation method to obtain the master equation for non Markovian 
interactions this is not important. 

We have discussed a number of derivations of the quantum jump approach so far. A 
different approach towards the description of single system dynamics has been proposed 
in (Teich et al, 1989; Teich and Mahler, 1992). In their method the dynamics described 
by the master equation is split into two distinct parts. One part changes smoothly the 
instantaneous basis of the density operator (coherent evolution) while the other part causes 
jumps between the basis states according to a rate equation. The instantaneous basis can be 
viewed as a generalisation of the dressed state basis. For a stationary state the basis states 
are fixed so that only jump processes occur. However, the approach is analytically quite 

43 



complicated for nonstationary processes and in addition there are interpretational problems 
(Wiseman and Milburn, 1993a). 

At this point we would like to explain briefly a recently proven connection (Brun, 
1996; Yu, 1996) between the quantum jump approach and a totally different concept, the 
Decoherent Histories formulation of quantum mechanics. A similar connection, although 
mathematically more involved, between the quantum state diffusion model and the Deco- 
herent Histories approach has also been established (Diosi et ai, 1995). The Decoherent 
Histories formulation of quantum mechanics was introduced by Griffiths, Omnes, and Gell- 
Mann and Hartle (Griffiths, 1984; Omnes, 1988; Omnes, 1989; Omnes, 1994; Gell-Mann 
and Hartle, 1990; Gell-Mann and Hartle, 1993). In this formalism, one describes a quantum 
system in terms of an exhaustive set of possible histories, which must obey a decoherence cri- 
terion which prevents them from interfering, so that these histories may be assigned classical 
probabilities. 

In ordinary nonrelativistic quantum mechanics, a set of histories for a system can be spec- 
iffed by choosing a sequence of times ti, . . . ,tiy and a complete set of projections {Paitj)} 
at each time tj, which represent different exclusive possibilities, i.e., they obey 



T.K,itj) = ^ , (113) 

K,{t,)Kdt,) = Sa^a'^a^it,) . (114) 

Note that the projection operators IP;^, are Heisenberg operators; one could represent them 
in the Schrodinger picture by 

K.itj) = e-^^'Pi/^K (115) 

The Schrodinger picture projection operators are assumed to be operators in the system 
space. 

A particular history is given by choosing one IP;^, at each point in time, specified by 
the sequence of indices {aj}, denoted a for short. The decoherence functional on a pair of 
histories a and a' is then given by 
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D[a, a'] = Tr^P^^M ■ ■ ■ Pl{h)p{to)\Pl>{h) ■ ■ ■ IP^^}, (116) 

where p(to) is the initial density matrix of the system. The decoherence criterion is now 
given by this decoherence functional D[h, h']. Two histories h and h' are said to decohere if 
they satisfy the relationship 

D[h,h']=p{h)6Hk', (117) 

where p{h) is the probability of history h. A set of histories {h} is said to be exhaustive 
and decoherent if all pairs of histories satisfy the criterion Eq. (117) and the probabilities 
of all the histories sum to 1. 

To establish a connection between quantum jumps and Decoherent Histories the idea 
is to use a system that interacts with the outside world in one direction. An example of 
such a system is a cavity. The counter outside the cavity is now modelled by a two level 
system that is strongly coupled to a bath so that both its coherence as well as its excitation 
is damped much faster than all time constants of the evolution of the system. One then 
defines the two projection operators 

IPo = i®|0)(0| IPi = i®|l)(l| . (118) 

These projections model the presence or absence of photons outside the system. It would 
be more general to consider more than one mode of the radiation field and the the proof can 
be generalised to that case. We now space these projections a short time 6t apart, and each 
history is composed of N projections representing a total time T = N6t. A single history 
is a string {ai, 0:2, . . . , o^n}, where aj = 0, 1 represents whether or not a photon has been 
emitted at time tj = {j — l)6t. Using this, it is possible to write the decoherence functional 
as 

D[/i,/i'] =Tr|lP,,e^^*(IP«,_,e^^*(...e^^*(IP«J^)(^|IP,;)...)IP,^} , (119) 

where C is the superoperator describing the time evolution according to the Bloch equations 
for the system (cavity) coupled to the two-level system. It is now possible to show that the 
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decoherence functional in fact obeys Eq. (117) to a very good approximation. It should be 
noted that the construction of the decoherent histories using the two operators in Eq. (118) 
closely resembles the derivation of the quantum jump approach as given by Hegerfeldt and 
Wilser (1991) and Wilser (1991). 

The crucial point in the quantum jump approach is the fact that we assume that we 
perform time resolved measurements on the photons that are emitted by the atom. These 
photons may be mixed with a local oscillator in a heterodyne detection as we did for the 
derivation of QSD but even there we assume time resolved measurements. A nice feature 
of the quantum jump approach has been that it allows us to describe the radiating system 
by a wavefunction instead of a density matrix. However, one may ask the question whether 
this is the only possible way to reach a wavefunction description of radiating systems. In 
fact it turns out that in some sense there is a complementary way to the quantum jump 
approach that also yields a wavefunction description. This method was proposed in (Hol- 
land and Cooper, 1996) and uses frequency resolved measurements instead of time resolved 
measurements. It turns out that again it is possible to decompose the density operator into 
pure states (Mollow, 1975) which, however, are now characterized by the number of photons 
that have been detected and for which the frequency instead of their precise emission time 
is known. 

C. Simulation of single trajectories 

After we have introduced and discussed different derivations of the quantum jump ap- 
proach we will now briefly explain how the quantum jump approach is used to simulate single 
quantum systems. We will describe the simulation approach for a decaying undriven cavity. 
The generalization to an arbitrary system should then be obvious. Carmichael (1993a) has 
given a precise relationship between the conditioned density operator contingent on a precise 
sequence of detection events (a "record" ) and the ensemble averaged density operator. He 
shows , for example, that if the zero temperature boson damping master equation is written 
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in Liouvillian form 

and we split the Liouvillian action £ as a sum of two terms, an anticommutator and a 
"sandwich" term, 



7 
2 

iC-S)+S , (121) 



Lp = a'a , p + ^yapa' 



then we may identify the "sandwich" term 5 as a jump operator. Equation (120) can be 
integrated formally as 

p{t) = exp{[{C-S)+S]t}p{0) 



v°°. /■* r*™ 


rt2 


> / dtm dt„,^i .. 


/ dt, 


m=0^0 ^0 


Jo 



where the quantity in curly brackets in Eq. (122) is labelled p^(t) by Carmichael and is the 
conditioned density operator describing a specific "trajectory" or detection sequence. We 
can write p^(t) in terms of the conditioned pure state projectors 

rcit) = mt)){'^cit)\. (123) 

The component exp[{C — S)At] propagates 'p^{t) for a time At without a decay being 
recorded: for the conditioned state |^c(^)) 

|W:(t + At)) = exp[-tH,sAt/h] |W:(t)) , (124) 

where the non-Hermitian effective Hamiltonian 

H^ff = H - ih-d^d (125) 

derives from the anticommutator in Eq. (121). Once a decay is registered, the gain in 
information is responsible for the jump 

47 



|^,(t))-^a|^e(t)). (126) 

The procedure adopted in quantum jump simulations can then be summarised as follows 
(Dalibard et al, 1992; Dum, ZoUer and Ritsch, 1992): 

1. Determine the current probability of an emission: 

AP = 7At(^|a^a|^). (127) 

2. Obtain a random number r between zero and one, compare with AP and decide on 
emission as follows: 

3. Emit if r < AP, so that the system jumps to the renormalised form: 

1^)— , \ • (128) 



4. Or no emission if r > AP, so the system evolves under the influence of the non- 
Hermit ian form 

^ — > ^ ^^^ ^^V -^^ ■ 129 

' ' (l-AP)i ^ ' 

5. Repeat to obtain an individual trajectory, or history. 

6. Average observables over many such trajectories. 

To reassure ourselves that this is all true, we note the history splits into two alternatives in 
a time At : 

(l^emit) with probability AP 

(130) 
\^noemissi<m) with probability 1 - AP 

Then in terms of the density matrix, the evolution for a step At becomes a sum of the two 
possible outcomes. 
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|^)(^| ^ AP|^emit)(^cmit| 

+ (1 - AP) l^no emit)(^no emit | (131) 

+ {1 - - HAt - ^ Ata^a}m(^\{l + - HAt - ^ Ata^d} 

~ |^)(^|-lAt[i7, |*)(*|] 

+ -At{2d\^){^\d^ -d^d\^){^\-\^){^\d^d} (132) 



so that 



— — = — T- [H, p] + — {2 apa^ — a^ap — pa) a} (133) 

as in the original master equation Eq. (120). We have now seen how the quantum jump 
approach can be used to simulate a master equation. In Section V we will see some examples 
of such simulations and also of single realizations of quantum trajectories. After this simple 
approach to simulating the master equation we now give a brief exposition of the idea 
of higher order unravellings of the master-equation (Steinbach et al, 1995a). To see the 
motivation we have to realize that the quantum jump approach is based on the simulation 
of the conditioned evolution of either a density operator or a state vector. However, at one 
point it is not a rigorous implementation of the trajectory concept. Because this method 
discretises time into small steps 6t , a quantum jump in the simulation takes a finite time 
6t , whereas in a simulation of quantum trajectories the information gained from detection 
should instantaneously be used in conditioning the quantum state of the system. This 
pinpoints the subtle difference between conditioned trajectories and the slightly simpler 
idea of evaluating the probability of decay quanta at discrete timesteps. The simplest way 
to remedy the fact that conditioning takes time in the simulation is to add evolution with 
the effective Hamiltonian to the projection step that has to be performed when a photon is 
detected. Having said this, the question arises at what point during the time interval 6t we 
need to condition the quantum state according to the result of the detection process. First, 
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it is worth noting that wherever we decide to do this, it would not change the accuracy of 
integrating the master equation in first order. Second, we may try to increase the accuracy 
by choosing a specific point in the interval 5t. Let us integrate the master equation to second 
order in 5t : 

psit + 5t) = psit) + ^ 5t ( [Cps\t + [Cps]t+6t ) + 0{5t^) . (134) 

Here C is the Liouville operator describing the evolution of the complete master equation. 
The terms that result from evaluating the right hand side of this equation can be cast into 
the following form (for details see (Steinbach et ai, 1995a)). By C we denote the reset 
operator that has to be applied after the detection of a photon. 

Ps{t + 6t)= Ups{t)U^ 

+ \5t UCpsit)C^W 
+ i(5t CUpsit)WC^ 

+ \5t^ UCCps{t)C^C^U^ + 0{5t^). 

(135) 

Here U denotes evolution under the influence of the effective Hamiltonian 

[/ = exp(-^iJeff5t), (136) 

which we call the "no-jump" evolution. The four terms on the right hand side of Eq. (135) 
represent four specific conditioned evolutions or mini-trajectories that the system might 
follow. An expansion into mini-trajectories is important because only then can the density 
matrix evolution Eq. (135) be simulated with pure states. The first mini-trajectory in Eq. 
(135) represents evolution without any jump, the second and the third represent a jump 
followed by evolution without jumps and vice versa respectively and the fourth includes two 
successive jumps followed by no-jump evolution. 
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We see that it is not sufficient to specify one point in the interval 5t at which to condition the 
density operator due to the quantum jump. We have to consider two points, at the beginning 
and at the end of 5t, and also the possibility of two immediately successive quantum jumps 
in order to increase the accuracy in 5t by one order. 

One can pursue this idea to obtain results which are accurate up to fourth order (in 5t). 
The master equation has to be integrated along the lines of a fourth-order Runge-Kutta 
method for ordinary differential equations. The result in fourth order then contains thirteen 
mini-trajectories (including the no-jump evolution) as follows 



ps{t + 6t) 



Uipsmi 

+ l5t UiCps{t)C^Ul + \5t CUipsit)UlC^ 

+ |5t UiCU2ps{t)UlC^Ul + |(5t U2CUips{t)U\C^Ul 

+ l6t^ Ui_CUiCps{t)C^UlC^U\ 
+ |(5t^ CUiCUips{t)UlC^UlC^ 



+ \5t^ Ui_CCUipsit)UlC^C^Ul 

2 2 2 2 

+ j^5t^ UiCCCps(t)C^C^C^Ul -\ 
+ j-^5e CCUiCps{t)C^UlC^C^ -\ 



^6t^ CUiCCps{t)C^C^UlC^ 
j^6t^ CCCUipsit)UlC^C^C^ 



+ j^Sf^ UiCCCCps{t)C^C^C^C^Ul + 0{6t^). 

(137) 

The subscripts on the non-Hermitian evolution U indicate the fraction of the time interval 
6t for which each particular U evolves the density operator, e.g., f/1/3 = exp (— -i Hcs St/3h). 
The way in which Eqs. (135) and (137) are turned into a Monte-Carlo simulation is clear: 
each mini-trajectory defines the conditioned evolution of the system and is assigned a spe- 
cific probability with which it occurs, analogous to the jump and no-jump probabilities in 
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the standard method. Just as in the standard procedure, a random number uniformly dis- 
tributed between and 1 is drawn to choose at random which of the mini-trajectories will 
govern the system evolution in the next timestep 5t. The no-jump evolution is tested first 
as this, for small 5t, is the most likely mini-trajectory. We note that the probability for 
evolution without detection remains unchanged as compared with the standard method and 
because the no-jump evolution is most likely the diversity of the mini-trajectories hardly 
influences the necessary computing time. However, if the no-jump mini-trajectory is not 
selected then one of the alternative trajectories in Eq. (137) must be chosen. For example, 
if the normalized state of the system at time t is |^(t)) then the state of the system after 
evolution corresponding to the fourth mini-trajectory in Eq. (137) is, 

which includes a renormalization factor A/". The state |\E') is evolved with the effective Hamil- 
tonian over two thirds of the timestep 6t using a fourth-order Runge-Kutta integration step 
(M0lmer et ai, 1993). After projection with the Lindblad operator C the evolution is con- 
tinued with the non-Hermitian Hamiltonian for the remaining third of the timestep 6t. Only 
then is the resulting state vector renormalized. The probability for this mini-trajectory to 
occur is given by the product of the factor 36t/8 and the renormalization A/". 

To illustrate the improvement the method of higher order unravellings presents (Stein- 
bach et ai, 1995a) one can simulate a laser driven two- level system using the ordinary 
first-order (fourth order integration of the effective time evolution operator) and compare 
the result to the same simulation using the method of higher order unravellings. In Fig. 9 
the simulation results of the inversion (0-3) = (|1)(1| — |0)(0|)/2 of the two-level system are 
plotted after 250000 runs and for a Rabi frequency Q = A equal to the Einstein coefficient of 
the transition, zero detuning and a time step 6t = 0.1 A~^. We clearly see that the first-order 
quantum jump approach (dashed line) deviates from the exact result (solid line) while the 
dotted line obtained from a fourth-order unravelling is much closer to the exact result. 
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D. A quantum system driven by another quantum system 

The next problem we want to investigate is that of a quantum system B driven by the 
radiation emitted from another quantum system A (Gardiner, 1993; Carmichael, 1993b). In 
the following we closely follow (Carmichael, 1993b). One could try to solve the problem by 
determining the dynamics of the driving system A first and from that the statistics of the 
emitted light. However, in general an infinite number of correlation functions is required 
to characterise the state of the light emitted from A. In semiclassical theory one could 
instead simulate the properties of the light by implementing a suitable stochastic process, 
unfortunately this is not possible in the quantum case. Therefore it is better not to divide 
the problem into two but to determine the dynamics of the composite system A (B B. To 
obtain the broken time symmetry one uses an interaction between A and B that is mediated 
by a reservoir R and one assumes the Born-Markoff approximation. A simplified version 
of the problem is illustrated in Fig. 8 where it is assumed that only one mode of each 
cavity needs to be assumed. The cavities have three perfectly refiecting mirrors and one 
with transmission coefficient T ^ 1. The Hamiltonians Ha and Hb describe the free cavity 
modes and any interactions that take place inside the cavities. Hji is the free Hamiltonian 
of a travelling wave reservoir R which couples the cavities in one direction only. The fields 
£^(0) and S{1) that couple to the cavities are written in photon fiux units. The complete 
Hamiltonian for A(B B (B R is 

H = Ha + Hb + Hr + Har + Hbr (139) 

with 

Har = ih {2ka)^^^ [aS\0) - S{0)a^] , 

HBR = lh{2KBY^^[bS^{l)-S{l)b^] , (140) 

where ka and kb are the cavity linewidths, and a and b are annihilation operators for the 
cavity modes. H describes two systems interacting with the same reservoir. It should be 
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noted that A and B couple to that reservoir at different positions in space. Usually spatially 
separated reservoir fields are treated as independent an assumption that cannot be valid for 
the geometry shown in Fig. 8 where the output from cavity A appears a time r = l/c later 
at the input of cavity B. The spatial separation of the two cavities can in fact be eliminated 
using the Born-Markoff approximation in the Heisenberg picture to relate ^^(0) and ^^(/). 
One obtains 

UAir)S^m\ir) = ^t(O) + i i2KAf' a , (141) 

where 

f/^(r) = e^(^^+-f^«+-f^^«)"/^ . (142) 

If x(t) is the density operator of system A © i? © i? we may define the retarded density 
operator of the system A © i? © i? as 

X'it) = UAir)xmiir) (143) 

and easily shows that x'{t) satisfies the Liouville equation with the Hamiltonian 

H' = Hs + Hr + Hsr, (144) 

where 

Hs = HA + HB + ih{KAtiBf^{a^b-ab^) , 
HsR = ih{[{2KAY/^a + {2KBy^^b\s\0)-H.c.} . (145) 

Now a and B couple to the reservoir at the same spatial location and in addition they also 

1 /9 

couple directly with coupling constant {kak^b) ■ Now one can easily derive the master 
equation for the density operator oi A® B p' = tr^jx'} and obtains 

p' = 1 [Hs, p'] + Cp'C^ - \c^Cp' - \p'C^C, (146) 

with 
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C = {2ka)^^^ a + {2kb)^^^ b . (147) 

Having found the master equation one can easily find the conditional time evolution that 
then allows to unravel the dynamics of the composite system A(B B. The time evolution of 
the conditional wave function lipdt)) between photodetections is governed by the Hamilton 
operator 



H = HA + HB-ih 



K,Aa^a + i^Bb''b + 2{K,Ai^BY^'^ab^] ■ (148) 



After a photodetection we have to reset the system using the operator C , i.e., I'lpd't)) ~^ 
C\ipc{t))- Now we are in a position to simulate individual trajectories of two coupled 
quantum systems. For applications of the theory , for example, to an atom driven by 
squeezed light or by antibunched light emitted from another atom see (Carmichael, 1993b; 
Gardiner, 1993; Kochan and Carmichael, 1994; Gardiner and Parkins, 1994). An early 
example of an investigation of atoms driven by antibunched light was discussed by Knight 
and Pegg (1982). 

E. Spectral information and correlation functions 

Until now we have discussed the quantum jump approach only in connection with quan- 
tities of the system or its resonance fluorescence that require a temporal resolution; no 
frequency information has been obtained as only broadband photon counting has been as- 
sumed. However, it would be nice to be able to use the quantum jump approach also for 
spectral properties of the resonance fluorescence of the system. Of special interest are , 
for example, the power spectrum of resonance fluorescence or the absorption spectrum of 
a weak probe laser. It turns out, that it is in fact possible to use the quantum jump ap- 
proach (Gardiner et al, 1992; Dum et al, 1992; M0lmer et al, 1993; Mu, 1994; Plenio, 
1994; Plenio, 1996; Hegerfeldt and Plenio, 1996) and also quantum state diffusion (Gisin, 
1993; Sondermann, 1995a; Brun and Gisin, 1996; Schack et al, 1996a), to calculate those 
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spectra. There are several different ways to derive quantum jump equations that enable us 

to calculate spectral information and in the following we will discuss three of them. 

One possible approach can be made via the method of quantum stochastic differential 

equations (QSDE) (Gardiner et ai, 1992; Dum et al, 1992). We already outlined the spirit 

of their approach in Eqs. (78) - (91). Following essentially Gardiner et al (1992) but using 

the notation used in Eqs. (78) - (91) one defines an output mode operator 

/•t+dt ( e^uj W^^l'^^ 

dM{t) = J^ ds-D2iY.'hi^^\ ekA^siakAe-'^"--"^^)^ (149) 

and then the spectrum in the Schrodinger picture as 

S(.) = lim imrK-,t)r(.,tmt)) _ ,150, 

t^°o t — Iq 

where we defined 

r{uj,t)= f'dM{s)e-"^'^'-'K (151) 

Jto 

We can now introduce the auxiliary wave function 

\P{t))=rico,tMt)) . (152) 

For these two wavefunctions one then obtains a Stratonovitch stochastic differential equation 
which can then be transformed into the Ito form. One then obtains 



d 



I \A.(+\\ \ ( _n u ..If, n \ / U(t)) I 



m)) 



iH^ff/n 



V m)) J 



(153) 



The spectrum is now by averaging 



s(.) . ii,„ mm (154) 

over many realization. Note that Eq. (153) for |0(t)) is just the conditional time evolution 
when no photons have been emitted. This equation is used to determine the jump times. 
The equation for \(3{t)) has a free evolution similar to \4>{t)) except for an additional rotation 
with frequency u. In addition it is driven by |0(t)). This driving can be interpreted as a 
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process where a photon is emitted into the modes described by dM{t) and where therefore 
the atom is deexcited. The approach given here can be used to simulate spectra in a number 
of situations (Dum, Zoller and Ritsch, 1992; Marte et al, 1993a; Marte et ai, 1993b). This 
approach is in fact similar in spirit to the approach of Schack et al (1996) who also assume 
the interaction of the system with one mode of the quantized radiation field. Schack et al 
then derive the quantum state diffusion equations that allow the calculation of the spectrum. 
It should be noted that using Eq. (153) and then following our derivation of quantum state 
diffusion from quantum jumps would have led to the same result. 

A slightly different approach to the problem more in the spirit of the derivation of the 
quantum jump approach of (Hegerfeldt and Wilser, 1991; Wilser, 1991) was undertaken in 
(Plenio, 1994; Hegerfeldt and Plenio, 1996). Their starting point is the photon number 
operator and they define the spectrum via the number of photons that have been emitted 
into certain narrow frequency intervals (Agarwal, 1974) or to be more precise, the number 
of photons in a certain mode of the quantized radiation field (an approach using the electric 
field operator following similar lines is also possible and is essentially equivalent to the 
approach of (Gardiner et al., 1992)). To remain close to a physical picture, they envisaged 
an experimental situation as depicted in Fig. 10. A part of the quantized radiation field (in 
a solid angle Qb) is observed by a broadband counter while in the rest of the space (solid 
angle Qs) a spectrometer (for example a Fabry-Perot) is situated. 

The broadband counter performs time resolved observations on the quantized radiation 
field. It is again assumed that the time resolved measurements can be modelled by a sequence 
of gedanken measurements that are performed in rapid succession as in the derivation of the 
quantum jump approach (Hegerfeldt and Wilser, 1991; Wilser, 1991). At a large time T, 
we then perform a spectrally resolved measurement of those photons that have entered the 
spectrometer. In that way one has measured a spectrum with a spectral resolution ~ 1/T, 
which is conditioned on the particular detection sequence one has found in the broadband 
counter. The derivation of the relevant set of differential equations to calculate the spectrum 
follows similar ideas as they were developed in the work of Hegerfeldt and Wilser (1991) 
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and Wilser (1991). The aim is to calculate differential equations for the conditional photon 
number in the mode kA which is given by 

tr{al^{0)a^x{0)p{t\h, ta, . . . , t„)} , (155) 

where u = |k|/c, k = k/|k|, and we assume that k is a vector pointing into the solid 
angle Qs of the spectrometer. The conditional spectrum of resonance fluorescence is now 
obtained by summing over all vectors k in Qs- For the calculation of Eq. (51) we need 
the state p{t\ti, t2, ■ ■ ■ , tn) where photons were detected at times ti, . . . , t„ only. With IPoq^ 
the projection operator onto the vacuum state for all modes with k G VLb and with the 
abbreviation 

A = \Pon,U{t, s^J n \Pm,U{sl, sl_,) , (156) 

fc=i 

where s^ are the times of those measurements where no photon was found, we flnd, for 

tn < t < tn+l, 

p{t\h,t2, ...)= Apitn + 0|ti, . . . , tn-l)A^ , (157) 

where p(t„ + 0|ti, . . . ,t„-i) is the state right after the detection of a photon in VLb', it is 
recursively given by 

p(t„ + 0|ti,...t„„i) = 

tm, {(1- iPof,,)f/(t„,sr„ijp(sr„ijti,...,t„_i)f/t(t„,c;iJ(i- iPof^j} ® iPon, . 

(158) 

Here trf7^{.} denotes the partial trace over all modes with a k vector that points into the 
solid angle VLb- As in the usual quantum jump approach, at this point the assumption enters 
that the photons detected in the broadband detector are absorbed during the measurement 
as in a real counter. One can show, however, that this assumption is not necessary for 
obtaining the equations of motion used in this section (Plenio, 1994). In fact one can just as 
well assume ideal quantum mechanical measurements instead and obtain the same results 
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(Plenio, 1994). This is of course a manifestation of the intuitive physical idea that photons 
emitted by the atom will not be reabsorbed as long as there are no reflecting mirrors close 
to the atom. Mathematically it is essentially a consequence of the Markov approximation 
which is in the approach of Gardiner et al (1992) incorporated elegantly in the formalism. 
The general principle is to find first the Heisenberg-Langevin equations for operators of the 
form 

Q\t) = UJ{t,0)\PonsQ\PonnUi{t,0) , (159) 

where Ui{t,0) is the interaction picture time evolution operator. With the abbreviation 

nin n-1 f nik ~j 

Vn ■■= n \PonAs]) n ^(^0+') n \PonAs^) , (160) 

j=0 fc=0 I 1=0 ) 

where Pong{t) and C{t) are defined similar to Eq. (75) without IPoqs, we then find for Eq. 

(55) 

tr{aL(0)o-(0)ajv,(0)p(t|ti, t2, ..•)} = (^I(aL^aM)° W^n) • (161) 

Having found the Heisenberg-Langevin equations the next step is to show that the occurring 
Langevin noise operators give a negligible contribution. Having shown this one arrives at a 
set of differential equations. As long as the solid angle covered by the broadband counter 
is not equal to 47r these equation still describe mixtures as the photons emitted into the 
spectrometer are not observed in a time resolved way. For the explicit form of the equations 
of motion see (Plenio, 1994; Hegerfeldt and Plenio, 1996). Going over to the limit of a 4 7r 
broadband counter yields equations that maps pure states into pure states and which are 
actually the same as those derived by Gardiner et al (Gardiner et al, 1992). The main 
ingredients in this derivation are the fact that the initial state of the quantized radiation 
field is the vacuum state and the Markov approximation. 

Using this approach it is straightforward to derive the quantum jump approach for the 
calculation of the absorption spectrum of a weak probe beam by an atom. The idea is to 
assume that the probe beam consists of a mode which is in a coherent state. 
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\a^x) = D{a^x)\0) , (162) 

where 

One then performs a unitary transformation that maps this coherent state onto the vacuum 
state. This leads to additional oscillating fields in the Hamilton operator and to an initial 
state in the probe mode which is now the vacuum. Therefore the approach outlined above 
may be used again. Without further complicated calculations one obtains the spectrum 
as the change in the photon number in the probe mode. It is possible to show that this 
definition of the absorption spectrum reduces to the stationary absorption spectrum for suf- 
ficiently long times (Plenio, 1994; Plenio, 1996). 

The approaches exhibited so far enable us to simulate spectral information (e.g. the spec- 
trum of resonance fluorescence). They were derived for both a definition of the spectrum 
via the electric field operator (Gardiner et ai, 1992) and via the photon number operator 
(Plenio, 1994; Hegerfeldt and Plenio, 1996). The physically motivated derivation of the for- 
malism by Hegerfeldt and Plenio (Hegerfeldt and Plenio, 1995b; Plenio, 1994) with a finite 
size spectrometer yields as a byproduct equations of motion for a system that is observed by 
a counter that does not cover the whole solid angle and/or has below unit efficiency. These 
equations of motion will later be used to illustrate the connection between the next photon 
and the any photon probability and to show that an inefficient photon counter will not 
measure the next photon probability but the any photon probability which is proportional 
to the intensity correlation function g'^'^^t). 

The approaches discussed so far gave equations that enabled us to calculate the spec- 
trum of resonance fluorescence directly from the norm of a component of a propagated 
'wave 'function. However, there is a third way to obtain the spectrum of resonance fluores- 
cence which is via the simulation of the correlation function (o'io(^ + ''")o"oi(t)) (where o"oi(t) 
is the Heisenberg operator corresponding to |0)(1|) and subsequent Fourier transformation 
of the simulation results. Indeed this is the way in which Dalibard, Castin and M0lmer 
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(M0lmer et ai, 1993; M0lmer, 1994) and others (Garraway et ai, 1995; Mu, 1994) employed 
the quantum jump approach to obtain simulations of the spectrum of resonance fluorescence. 
The simulation procedure for a correlation function of the form C(t + t) = {A{t + r)B{t)) 
runs as follows. First one evolves using the MCWF approach a wave function |0(O)) towards 
|0(t)). Then one forms the auxiliary wave- functions 

|x±(O)) = -^(l±i?)|0(t)) , (164) 



|x'±(O)) = ^(l±^i?)|0(t)) , (165) 

where the /i±,/i'^ are normalization constants. Now one has to evolve each of these four 
wave functions according to the MCWF procedure and then to form 

c±{r) = {x±{r)\A\x±ir)) , (166) 

4(r) = (x±(r)|A|x'±(r)) , (167) 

from which on obtains 

C{t,T) = - fj,+c+{r) — ji^c^ir) — in'^c^ij) + ifi'_c_{T) . (168) 

It can be shown that this procedure produces the correct ensemble averages (Dalibard et 
ai, 1992). A subsequent Fourier transform of the simulation results yields a spectrum. This 
procedure has been used , for example, (Garraway et ai, 1995) to simulate the spectrum of 
a three level system in a V configuration. However care has to be taken in this procedure 
because if one only performs a single Fourier transform to obtain the stationary spectrum 

Si{co):=2Re rfre^^"((7io(r)(7oi(0)),, , (169) 

Jo 

one inevitably runs into difficulties as one can only simulate finite times. Then, however, the 
Fourier transformation yields spurious negativities in in the power spectrum which are due to 
finite time effects. This problem can be circumvented by using a finite time double integral 
of the correlation function as required in the definition of the time dependent spectrum. 
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St{uj) := i f^dt' r dt"e'^^''-'"\a^o{t>oi{t")) ■ (170) 

1 Jo Jo 

This quantity has the advantage that it is manifestly positive for each reahzation as it is of 
the form (A'^A). However this considerably slows down the simulation procedure because 
the integral transform is now much more complicated. 

A different point of criticism should be mentioned concerning the interpretation of a 
single run of this simulation procedure. While the derivation of the simulation procedure 
by Hegerfeldt and Plenio (Hegerfeldt and Plenio, 1996) clearly shows that a physical in- 
terpretation of a single run in the schemes derived in Hegerfeldt and Plenio (1996), Plenio 
(1994) and Gardiner et al (1992) is possible, although this is not obviously the case for the 
simulation procedure of M0lmer et al (1993) . A problem lies in the fact that four different 
wave functions have to be propagated in the procedure, each of which might follow different 
jump sequences. Furthermore even if only one jump sequence is followed, it is not clear 
what the interpretation of the correlation function is. This may be illustrated by the fol- 
lowing example. We now would like to simulate a trivial correlation function of the form 
(l(t + T)A(t)). We would expect that in each realization the simulation procedure would 
give the same result as for the single time expectation value {A(t)). This, however, is not 
the case as is easily seen in the following example. We assume a spontaneously decaying 
TLS and the initial state 

\m) = {-^)'^\\0) + e\l)) . (171) 

The emission-free time evolution is given by 

f/o(t + t',t) = e^*' , (172) 

where 



M 



^0 ^ 

y —T22 I 



(173) 



If we assume that no jump at all occurs until t + r we obtain according to the procedure of 
M0lmer et al (1993) the single run result for the correlation function (l(t + r)o'oi(t)) to be 
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(l(t + r)aoi(t)) = ^^^L_ (174) 

while for the expectation value of the operator o'oi(t) at time t = l/Fn is given by we obtain 

{aodt = l/Tn)) = l . (175) 

These two expressions obviously differ. Therefore we can conclude that the simulation proce- 
dure of M0lmer et al (1993) yields the correct ensemble results; however if we are interested 
in questions concerning single runs or spectra conditioned on a given jump statistics care has 
to be taken and it is safer to resort to the approaches in (Gardiner et al., 1992; Dum et al., 
1992) and Hegerfeldt and Plenio (Hegerfeldt and Plenio, 1996; Plenio, 1994). Applications 
to conditioned spectra will be given later in Section V. 

Here we have only described in detail quantum jump approaches to calculate spectral 
information. Of course one can also obtain a quantum state diffusion simulation of spectral 
information. One can derive the equations either by starting out from the quantum jump 
equations given above and then follow our derivation of the quantum state diffusion equations 
as a result of heterodyne detection. Otherwise one can deduce the equations by assuming 
that the system couples to an additional output mode (i.e. similar to (Gardiner et al., 
1992) but here literally one mode of the quantized radiation field is chosen) and then writes 
down the master equation for that enlarged system. From that one can then easily obtain 
quantum state diffusion equations (Brun and Gisin, 1996). However, if one tries to apply the 
equations to the calculation of time correlation functions of the form {A{t+T)B{t)) one again 
encounters interpretational problems similar to those explained for the approach of M0lmer 
described above. It should be noted that the approach to calculate time correlation functions 
given by Gisin (1993) does not lead to {A{t + T)B{t)) as usually defined in quantum optics 
but to a different correlation quantity which is difficult to interpret (Sondermann, 1995a). 
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V. APPLICATIONS OF THE QUANTUM JUMP APPROACH 

A. Photon statistics 

In the previous sections we have introduced and discussed the quantum jump approach 
illuminating many different approaches to it. After these sometimes formal considerations we 
would like to give a number of examples to give a better understanding of the formalism and 
its physical implications. The examples will be drawn mainly from two physical situations; 
single trapped ions driven by lasers and electromagnetic fields in cavities, i.e., cavity QED. 

We start by illustrating the difference between single realizations and ensemble aver- 
ages by investigating a laser driven two-level system (TLS) damped by a zero temperature 
quantized radiation field. The master equation is then given by 

p = --[HA,p\-V{aiip + paii] + 2TaQipaiQ , (176) 

Ti 

with 

fiQ. 

i7^ = -;^A|l)(l| + — (|0)(1| + |1)(0|) (177) 

where 2T equals the Einstein coefficient of the observed TLS, Vt the Rabi frequency and A 
the detuning of the laser. 

Assuming standard broadband photodetection, the quantum jump approach gives for 
the conditional time evolution between detections 

m)) = e-^=//*/^|^(0)) , (178) 

with 

Heff = HA-thT (179) 

and for the normalized state after the detection of a photon 

mu)} = '^^ = 10) . (180) 
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The probability for not having a jump in the time interval [0,t] if the initial state is the 
ground state, \ip{0)) = |1) is given by 

A2 — Al A2 — Al 

with 



Using the parameters fi = 5 F and A = the upper state population pu of the ensemble 
evolves as shown in Fig. 11 where the initial state is the ground state. The population 
rises from zero, undergoing some Rabi oscillations and then tends towards a steady state. 
So much to the ensemble average. Now let us look at individual realizations of the time 
evolution using the same parameters and initial conditions. The result for one possible 
realization is shown in Fig. 12. The picture is strikingly different from Fig. 11 in that 
the time evolution is not smooth but exhibits jumps and it does not tend towards a steady 
state. We rather observe that initially the system starts to perform a Rabi oscillation. As 
the population in the upper level grows in time so does the probability for an emission of 
a photon. This oscillation is then terminated by the emission of a photon which brings the 
atom back to its ground state. Then the whole process starts again. Figs. 11 and 12 show 
little similarity; however, averaging over many individual realizations shown in Fig. 12 lead 
to a closer and closer approximation of the ensemble average. To illustrate this we have 
plotted the ensemble result together with the average over A^ = 100 and 10000 individual 
realizations in Fig. 13. It is found that the root mean square deviation of the simulated 
average from the exact ensemble result are of the order of \/yN . We have examined the 
photon statistics of a two level system by detecting all the photons emitted by the atom 
with high time resolution as we assumed that the counter has unit efficiency and covers the 
whole solid angle. Assuming this, one would be able to measure the next photon probability 
density /i(t). However, in a real experiment the counter efficiency is much less than unity 
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due to imperfect quantum efficiency and because only a finite solid angle is covered by the 
counter. Therefore the question is which function we actually measure when we perform 
a real experiment in which we determine the detection time of the 'next' photon in our 
imperfect counter. 

To answer this question we would like to employ the quantum jump approach assuming a 
perfect photon counter which, however, does not cover all space. If it covers a solid angle Qb 
then this is equivalent to a An counter of efficiency f3 = ^2^/47?. The equation of motions for 
this setup up have been derived in section IV in connection with the simulation of spectral 
information (see Eqs. (155)-(161)). The efficiency (3 of the detection process is the fraction 
of the emitted photons that are actually detected. We obtain the conditional equation of 
motion for the density operator under the assumption that no photon has been detected in 
the counter, 

P = -Jl {HeffP - pHeff} + 2r (1 - P) aoipaio . (183) 

Heff is given by Eq. (179) and one should note that now the conditional time evolution does 
not map pure states onto pure states. This is of course due to the incomplete information 
gained from the imperfect counter. Therefore we have to average over all possible events 
that the counter could not detect and this leads to a mixture. One observes two familiar 
limits: for zero counter efficiency we recover the ordinary optical Bloch equations, whereas 
for unit efficiency f3 = 1 we find the effective time evolution inferred from the assumption 
that no photon has been found in the whole solid angle given a perfect counter. 

Having found Eq. (183) we can now calculate the probability that the detector finds no 
photon in the time interval [0, t] 

Po{t)=tr{po{t)} , (184) 

which reduces to Eqs. (56) or Eq. (67) in the limit (3-^1. The next count rate is then the 
negative time derivative of Eq. (184) 

/i,^(t) = 2r/?pn(t) , (185) 
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where P22 is obtained by solving Eq. (183). Instead of solving this equation analytically we 
plot the solution for several different counter efficiencies 13. We assume F = 1 and f] = 5 F. 
In Fig. 14 we plot the next detection probability for counter efficiency /3 = l,/3 = 0.1 and 
j3 = 0.0049 where the latter value of (3 is computed from Mandels antibunching experiments 
(Kimble et ai, 1977; Dagenais and Mandel, 1978). We observe that for decreasing counter 
efficiency the function Eq. (185) approximates the "any photon" probability of the whole 
ensemble, which becomes more and more closely proportional to the intensity correlation 
function g^^\t). So far we have we have shown how to calculate the next detection rate for 
different counter efficiencies from the quantum jump approach. We have plotted the results 
numerically. However it is also possible to derive an analytical expression in a somewhat 
different way. This was done in Kim et al (1987) and was in fact one earlier approach to 
the problem. As described earlier one can utilise the fact that there is a very simple relation 
between the "any photon" rate //3(t) of the complete ensemble and the next detection rate 
/i^^(t) for a counter with efficiency f3. This relation can be found using 

PiiW = Pu\t) + fdt'puit - f) 2TPp[^,\t') , (186) 

Jo 

and 

hAt) = 2^PPn\t) , (187) 

I^it) = 2TPpuit), (188) 

where Pu (t) is calculated from Eq. (183). Inserting this and taking the Laplace transform 
this yields 

P 1 - h,p{z) 
The intensity correlation function g^'^\t) is related to the any photon rate Ip{t) for an 
imperfect counting process by 

2Vp^^{^)Pg^''\t) = Ip{t). (190) 
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Therefore we obtain 

where g^'^\z) is the Laplace transform of the intensity correlation function g^'^\t). 

One should note that the intensity correlation function does not depend on the efficiency 
f3 as it is normalized with respect to the stationary detection rate. From expression Eq. 
(191) one observes that, up to a normalization factor, the next photon probability tends 
towards the intensity correlation function. Inserting the well known expressions for the 
stationary state of the TLS and the intensity correlation function for zero detuning 

PrA.co) = j^j3^ (192) 



and 



^ ^^^^ z{{z + 3T/2)^ + n^-Ty4) ^^^^^ 

into Eq. (191) we obtain 

^''^^'^ ^ {z + r){z{z + 2r) + Q^{z + f3r)) ' ^^^^^ 

This Laplace transform can be inverted easily. Plotting this function exactly reproduces 
Fig. 14 if one normalises the maximum of Ii^p{t) to unity. 

B. Intermittent fluorescence 

After this discussion of the photon statistics of the resonance fluorescence of a two-level 
system we now proceed to illustrate the formalism for a more complex system, namely, 
a three- level system in V-configuration as shown in Fig. 1. We assume that the <-> 
1 transition is strongly allowed while the ^^ 2 transition is metastable. In a typical 
experiment one would have a lifetime of the order of seconds for the metastable level 2, 
while the unstable level 1 has a lifetime of several nanoseconds. The system is irradiated by 
two lasers, one on each transition. The Rabi frequency Vti of the laser driving the ^^ 1 
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transition is assumed to be much larger than the Rabi frequency ^2 on the ^^ 2 transition. 
If one observes the intensity of resonance fluorescence on the strong transition under these 
conditions one typically obtains a result as shown in Fig. 2 or in a more schematical 
representation in Fig. 4. Long periods of brightness with many photon counts (bright 
periods) are interrupted by prolonged periods with no photo-detections (dark periods). As 
we discussed previously a simplified treatment of this situation using rate equations was 
given by Cook and Kimble (Cook and Kimble, 1985; Kimble et ai, 1986). However, as 
the described situation involves lasers, a more detailed treatment is required as coherences 
can play a crucial role in the time evolution of the system. Such treatments were initially 
undertaken using Bloch equations but this is not the most natural description of the problem. 
Such a natural description of the problem was provided by the quantum jump approach. 
In the following we will use the quantum jump approach to calculate the photon statistics 
of the V-system and we will use it to gain interesting and sometimes surprising insights 
into the single system dynamics. It will turn out that the rate equation treatment of Cook 
and Kimble, while giving a qualitative picture, is insufficient to explain many interesting 
properties of the system. A similar analysis can be carried out for other systems such 
as a A configuration (Agarwal et al, 1988b; Hegerfeldt and Plenio, 1995a; Plenio, 1994), 
and in a system where both upper levels couple strongly to the ground level but are close 
together (Hegerfeldt and Plenio, 1992; Hegerfeldt and Plenio, 1993; Hegerfeldt and Plenio, 
1994; K5hler, 1996). 

In the following analysis we will always assume that both bright and dark periods are 
much longer than the lifetime of the unstable level 1. This condition is necessary from a 
physical point of view as otherwise one would not be able to distinguish between a dark 
period and the time interval between two successive emissions in a bright period. Also a 
bright period consisting of approximately one photon has little meaning. With the detuning 
Aj and Rabi frequencies Qi of the lasers on the ^^ i transition we obtain the condition 
(Hegerfeldt and Plenio, 1995b) 
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^^ ^^ i rf, + (A, - A,y • ^^^^^ 

where 2rii equals the Einstein coefficient of level 1. This condition assumes that the Einstein 
coefficient of level 2 is negligible. Again this can be cast into an analytical form. We require 
(Hegerfeldt and Plenio, 1995b) 

i Zi i Zo-L 11 / X 

'' ^ 16Air?i + {fif + 4A2(Ai-A2)}2 ^^^^^ 

which expresses the fact that spontaneous emissions from level 2 are much less frequent than 
stimulated transitions. This also implies that there are sufficiently many long dark periods. 
To be able to calculate the mean lengths of bright and dark periods we need to define 
precisely what we mean by bright and dark periods. To distinguish between bright and 
dark periods we introduce a time Tq. If our perfect photo detector fails to detect a photon 
in a time interval [0, To] where it has found a photon at time t = we speak of a dark 
period of length T^. If, however, in a time interval [0,T^] the time between successive 
photon detections is always less than Tq we have a bright period of length T^. Using these 
definitions and the next photon probability density /i(t) we obtain for the mean length of a 
dark period 

while the mean length of a bright period is 



1 /o"Wt7i(tO 



tl{to) = ^^ ':zyrZ. ■ (19^ 



These expressions appear to be quite complicated but they can be simplified when we use 
the assumption that bright and dark periods are long compared to the lifetime of level 1. In 
that case it is obvious that the time evolution of the system has two widely different time 
scales, one determining the emission rate during bright periods, the other giving the rate at 
which dark periods occur. This is reflected in the fact that the probability to flnd no photon 
in the interval [0, t] has a slowly decaying tail as we can observe from Fig. 15 where we plot 
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Po{t) for the parameters Qi = 2rii,r22 = O.SSFii, Ai = A2 = and r22 = 0. If we chose 
To such that it is much larger than the mean time between photon detections in a bright 
period while still being much shorter than the mean length of a dark period we can reliably 
distinguish between bright and dark periods. The choice of Tq implies that 



Po{To)=pe 



-2To/mAi 



P< 1 



(199) 



where Ai is the smallest eigenvalue of the effective Hamiltonian of the system 



H, 



-f 


^2 
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-f ^Tn + Ai 
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A2 
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e// 



Simplifying Eqs. (197) - (198) we then obtain 

1 



(200) 



ToiTo) 



-2/mAi 



(201) 



and 



Tl(To) = ^ 
P 



(202) 



tl{To)= dt't'hit')/ dt'hit' 
Jo JO 



(203) 



where Tl is the mean time between successive photo-detections in a bright period 

dt't'hit')/ / 
/o Jo 

The mean time interval between successive emission in a bright period can be found from 

Eq. (203) or more easily from 

1 



Tl 



(204) 



^J- llPll 

where pf/"^ is the population of the upper level 1 of a TLS driven by a laser of Rabi frequency 
VLi = 2 Til and Ai = 0. Now it is quite easy to obtain analytical expressions for T^ and Tj,. 
We find 

16A2r2, + (l]2_4A2(A2-Ai))2 



2i '-|^i '2I 11 
"^ 2Vli + 2{Ai-/\2? "" 



(205) 
(206) 
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We see that under the conditions Eqs. (195-196) both bright and dark periods are much 
longer than the hfetime of the unstable <-^ 1 transition. 

Investigating the average rate ^ ^j, at which we observe quantum jumps, i.e., the onset 
of dark periods, we see that it depends on the detuning of the laser on the <-^ 2 transition. 
If we assume that the strong laser {Qi ^ Fi) is resonant (Ai = 0) we observe that the ratio 
^ ]_^ becomes minimal when A2 = and maximal for the Autler-Townes or Stark split 
detunings A2 = ±^. 

1 ^ mlQlT^, r^, + Aj 

Tn + n 16 AlHi + [nl - 4A1)2 AVl, + 2 Ai + f]f ' ^ "" ^ 

This dependence on the detuning A2 reflects the fact that due to the strong driving of the 
^^ 1 transition the lower level exhibits Autler-Townes splitting. The two effective levels 
are shifted by ±r2i/2 with respect to the original level 0. To obtain long dark periods 
one needs to bring the weak laser on the ^^ 2 transition into resonance with one of the 
dressed states of the ^^ 1 transition. The resulting frequency dependence is shown in 
Fig. 16. Before we proceed with the investigation of the single system behaviour, we show 
again, now quantitatively, that the existence of bright and dark periods in the resonance 
fluorescence has a visible effect on the ensemble quantities, too. To see this, we plot the 
intensity correlation function g^'^^r) in Fig. 17. We observe that there is a slowly decaying 
contribution for times r ^ F^2^ quite similar to the behaviour of the next photon probability 
density /i(t) shown in Fig. 15. One can show that in this regime we can approximate g^'^^r) 
by 

gi2) (r) = l + ^exp I (^^ + 1-) r } , (208) 

where To and T^ are given by Eqs. (205)- (206). Therefore from the measurement of the g^"^^- 
function we can infer the existence of long bright and dark periods and one can determine 
their mean lengths. The fact that this information is hidden in the intensity correlation 
function is not surprising as it is proportional to the a photon rate /(r) which is related to 
the next photon rate Ji (r) via 
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/(r) = /i(r)+/ rfr7(r-r')/i(r') . (209) 







Later we will also show that the existence of bright and dark periods also leaves its finger- 
prints in the spectrum of resonance fiuorescence as well as in the absorption spectrum. 

After this analytical discussion of the photon statistics of the V-system we will now 
investigate how the system time evolution for a single realization will typically look like. 
For the parameters Qi = 2rii,i72 = 0.2rii, Ai = 0, A2 = and r22 = we have plotted 
both the time evolution of the population pu of the unstable level 1 and population P22 
of the metastable level 2. In Fig. 18 one observes that in some regions a rapid change of 
Pii accompanied by many detections takes place. However, there are also regions where no 
photon is found and most population is found in level 2 while there is still a remnant in level 
1. One should note that, as seen in Fig. 19 the population in level 2 grows continuously and 
does not jump from level into level 2 as one might expect from a rate equation picture. 
Nevertheless there is a jump, but not at the beginning of the dark period, but at its end 
(Wilser, 1991; Garraway, Knight and Steinbach, 1995) as we can see in Fig. 20 where we 
have magnified the time evolution of pn in a dark period. The jump occurs after a long time 
and marks the end of the dark period as after the jump the population is in the ground state 
which is strongly coupled to state 1. It should be noted that in the simulation shown in Figs. 
18-20 the decay rate of the metastable level was assumed to be r22 = which implies that the 
dark period always ends with an emission from level 1. The continuous change of population 
towards the metastable state in a dark period nicely clarifies the importance of the failure to 
detect photons (null measurements) for the time evolution of the wavefunction. The failure 
to detect a photon provides us information about the system state which is described by the 
wave function. In this case the no-detection event tells us that we are more likely to find 
the system in the metastable state. Accordingly the wavefunction quickly (within several 
lifetimes of the unstable level) tends towards the metastable state if we have not found a 
photon for this time. The evolution is continuous as we assume (in the limit of the Markov 
approximation) continuous measurements of the radiation field. Therefore the failure to 
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detect a photon in an infinitesimal time interval [t, t + 5t] leads to an infinitesimal change 
proportional to 5t. On the contrary the detection of a photon leads to a discontinuous change 
of the wave function because one quantum of energy has leaked out of the system into the 
environment. Therefore the system has to change disco ntinuously and jumps back into its 
ground state. Again this change of the wavefunction is due to our increased knowledge about 
the system. 

C. Prom quantum jumps to quantum state diflFusion 

So far we have only illustrated the quantum jump approach to the photon statistics of a 
single ion. In Section IV we have seen that other pictures are possible, namely, the quantum 
state diffusion picture. Instead of presenting a large number of examples for the quantum 
state diffusion model we rather want to illustrate how the transition from the quantum 
jump picture to the quantum state diffusion picture takes place. In Section IV we have 
seen how to derive the quantum state diffusion equations from a quantum jump description 
of a balanced heterodyne detection experiment. In this derivation we had to assume the 
limit in which the photon number in the local oscillator tends to infinity. In the following 
we will illustrate this limit by choosing a number of finite values for the photon number of 
the local oscillator. To keep the following short analysis as simple as possible we in fact do 
not consider the case of balanced heterodyne detection but of homodyne detection (Vogel 
and Welsch, 1994). Instead of a cavity we will investigate the time evolution of a laser 
driven two-level system with upper level 1 and lower level 0. We follow the presentation 
give in (Granzow, 1996). For the simulations we need to know two quantities. The Lindblad 
operator for the homodyne detection scheme is given by 

L=|0)(l| + al (210) 

where a is the field amplitude of the local oscillator. As the Lindblad operator is changed 
and depends on a so does the effective Hamilton operator between photon detections. We 
have 
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iJe// = -^A|l)(l| - ^(|1)(0| + |0)(1|) - thVniL^L + 2a*L + \a\^) (211) 

where 2rii equals the Einstein coefficient of the upper level 1, Q is the Rabi frequency on 
the ^-> 1 transition and A the detuning. The probability for a jump in the time interval 
6t is given by 

Pc = 2Tn [{L^L + a*{L) + a{L^) + |ap] 6t (212) 

These expressions can be derived analogously to the procedure that we applied in the de- 
scription of the balanced heterodyne detection. In Fig. 21 we have plotted single realizations 
of a two-level system with Rabi frequency Q = AT and local oscillators amplitudes a) a = 0, 
h) a = 0.5, c) a = 1 and d) a = 10. One clearly observes that with increasing a the number 
of jumps increase while their amplitudes decrease. For a = 10 we already see a behaviour 
very close to that one would obtain in the limiting case a = oo. 

D. A decaying cavity 

So far we have illustrated the quantum jump approach in the context of single ion 
resonance fluorescence. Now we would like to discuss a different kind of problem, namely, 
cavity QED (Haroche, 1984). That means we are interested in the field states inside a 
cavity as well as the time evolution of an atom interacting with such a cavity (Grochmalicki 
and Lewenstein, 1989b; Imamoglu, 1993). Again these problems can be formulated in the 
framework of the quantum jump theory because both the losses of the cavity to the outside 
world as well as spontaneous emission of an atom give rise to a master equation of Lindblad 
form which can be simulated using wavefunctions. A derivation of the simulation equations is 
possible using the physical picture of photon counters surrounding the system. For example, 
a broad band photon counter outside the cavity will detect photon losses of the cavity. As 
an initial example to illustrate the physical insight that we gain from the quantum jump 
approach, we would like to investigate the time evolution of a field state in a lossy cavity. 
As an initial state we chose the state 
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w) ^ i^m . (213) 

In an appropriate interaction picture, the time evolution under the condition that no photon 
has been found by the photon counter outside the cavity is given by the effective Hamihonian 

H^ff = -iTa^a (2f4) 

while the normalized state after the detection of a photon is given by 

Ife) = *U . (215) 

'{ij\a^a\ij) 



In Fig. 22 the mean photon number for a single realisation is shown. Initially no jump takes 
place and the effect of this failure to detect a photon outside the cavity is that the mean 
photon number inside the cavity decreases; it becomes more and more unlikely to find a 
photon inside the cavity because if there was a photon in the cavity then it would leave the 
cavity, leading to a photon count. However, in the simulation shown in Fig. 22 we finally 
observe a photon outside the cavity and at that moment we know that there have been 
photons in the cavity. Calculating the state after the detection, given that we started with 
the superposition Eq. (213) we find that the post-detection state is a Fock state containing 
9 photons. This implies that the mean photon number after the photon detection is actually 
higher than before the detection. After this first detection, the cavity continues to decay and 
now each detection of a photon outside the cavity decreases the mean number of photons 
inside the cavity, while now the number of photons remains invariant under the conditional 
time evolution from the relevant Fock states. 

In the previous example we saw that the exponential decay of the field mode in the 
ensemble average is the result of the superposition of many single realizations in which the 
cavity excitation changes discontinuously at random times. Coherence between component 
parts of a superposition state changes in amplitude and is eventually destroyed. However, 
this is not the only way in which coherence between superposition states is destroyed. In 
the following example we show that the decay of coherence of a Schrodinger cat of the form 

76 



|^) = (|a) + |-a)))/||.|| (216) 

is due to a randomisation of the relative phase between the two coherent states while the 
modulus of the relative phase remains invariant under the no-jump evolution (Garraway and 
Knight, 1994a). In each individual realization the cat state remains a cat state. Although 
the amplitude of the two coherent states decays the relative phase between the two coherent 
states just changes its sign. In fact if the state before the detection is given by Eq. (216) 
then after the detection of a photon we have the state 

l^+) = «l^)/IMI 

= (!«)- I -«)))/IMI ■ (217) 

The conditional time evolution when no photon has been found is given by 

|^(t)) = (|ae-^(*-*')) - I - «e-^(*-*'))))/||.|| (218) 

so that the amplitudes of the coherent states decay exponentially while their relative phase 
remains unaffected. Averaging over many realizations, however, leads to a decaying relative 
phase as random phases tend to cancel out. It is interesting to note that the rate of decay 
for the relative phase of the Schrodinger cat is given by r|ap while the amplitude decay rate 
is given by T. Therefore in the ensemble average the cat decoheres before the amplitudes of 
its constituents are significantly affected (Garraway and Knight, 1994b). 

Another example that shows that the loss of coherence can be due to phase shifts at 
random times is the phenomenon of revivals in the Jaynes-Cummings model (Shore and 
Knight (1993) and references therein). In this phenomenon a two-level atom is initially in 
its excited state 1 while the cavity field is in a coherent state with amphtude a, 

1^) = \1){1\ (g)\a){a\ . (219) 

This initial product state becomes entangled (Ekert and Knight (1995) and references 
therein) and the average excitation of the atom rapidly tends towards 1/2. The reason 
for this is the fact that the frequency of Rabi oscillations depends on the number of photons 
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in the field mode (Narozhny et al, 1981). The different Rabi frequencies quickly decohere so 
that one observes an average excitation of 1/2 of the atom. However, as the frequencies are 
discrete they partially rephase after a longer time and the excitation of the atom rises again. 
This revival is shown , for example, in Fig. 23. Revivals of this kind have been observed 
experimentally by (Meekhof et al., 1996; Brune et al, 1996) and of the micromaser kind by 
Rempe et al (1987). If the cavity is damped, however, these revivals are much weaker and 
will vanish for strong damping as shown in Fig. 24 where we have chosen the same param- 
eters as in Fig. 23 but with an additional cavity damping constant of 7 = O.Ol^', where g is 
the atom-field coupling constant. We assumed that the atom does not decay spontaneously. 
No substantial revivals can be observed (Barnett and Knight, 1986). What will we see for 
an individual realization? If we assume that the atom does not decay spontaneously we 
obtain the effective Hamiltonian 

Heff = Hsys — ih'ya'a (220) 

that generates the conditional time evolution of the atom cavity system if no photons are 
detected outside the cavity. Hgys generates the free evolution of an undamped atom cavity 
field system while the second term on the right hand side describes the damping of the cavity 
at a rate 27. The state after the detection of a cavity photon which is 

l'^.) = f[ . (221) 

A single realization for the parameters of Fig. 24 is shown as in Fig. 25. Now the revivals 
persist, however, each jump that occurs at a random time introduces a phase shift. An 
average taken over many individual realizations then leads to a quick decay of the coherence 
and revivals are not observed any more. A similar analysis can be performed for a lossless 
cavity with an atom that may spontaneously decay (Burt and Gea-Banacloche, 1996). 

In these examples we only considered undriven atoms in a cavity. Driven two-state 
systems inside a cavity can of course also be investigated using the quantum jump approach 
(Alsing and Carmichael, 1991; Tian and Carmichael, 1992) but we do not discuss this in 
detail here. 
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E. Other applications of the quantum jump approach 

So far we have illustrated the quantum jump approach for a number of examples in 
single ion physics and cavity QED. The quantum jump approach can and has been applied 
to a large number of problems which we cannot discuss in detail here. Amongst these 
examples are the discussion of correlations between quantum jumps in two stored ions. Such 
correlations were observed in an experiment (Sauter, Blatt, Neuhauser, and Toschek, 1986) 
in which the intermittent fluorescence of two Ba^ ions in an ion trap were investigated and 
correlations were found between the quantum jumps in the two ions exceeding those expected 
from independent jumps. However, Itano et al (1988) could not confirm these correlations 
in their experiment (Itano et al, 1988). Theoretical investigations (Agarwal et al, 1988a; 
Hendriks and Nienhuis, 1988; Lawande et al, 1989a; Lawande et al., 1989b; Lewenstein and 
Javanainen, 1987; Lewenstein and Javanainen, 1988) predicted that correlations should be 
small if the ions are separated by many optical wavelengths. If the ions are closer than an 
optical wavelength together, correlations might be observable. For such small separations 
even two two-level ions could exhibit quantum jumps (Yamada and Berman, 1990; Kim et 
al., 1989; Lawande et al., 1990). 

Another example of the application of the Quantum jump approach is the investigation 
of the Quantum Zeno effect (Misra and Sudarshan, 1977). The Quantum Zeno effect was 
measured in an experiment by Itano et al (1990) which was originally proposed by Cook 
(1988). The experiment used a two-level system. The ground state population was measured 
via coupling the lower level strongly to a rapidly decaying third level. Observation of reso- 
nance fluorescence then indicates that the system was in its ground state. This experiment 
was investigated theoretically independently by Power and Knight (Power, 1995b; Power, 
1995a; Power and Knight, 1996), by Beige and Hegerfeldt (Beige and Hegerfeldt, 1996a; Beige 
and Hegerfeldt, 1996b) and by Mahler and coworkers (Mahler and Weberruf5, 1995) using the 
quantum jump approach and, e.g., by Frerichs and Schenzle (1991) using Bloch equations. 
The investigations using the quantum jump approach helped to understand the experiment 
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from the single particle point of view and will be important in the analysis of future Quan- 
tum Zeno experiments (Plenio et ai, 1996) using single ions instead of around 5000 in the 
experiment by Itano et al (1990). 

Cohen- Tannoudji and coworkers have applied the quantum jump approach to the problem 
of lasing without inversion (Cohen- Tannoudji et ai, 1992a; Cohen- Tannoudji et ai, 1992b; 
Cohen- Tannoudji et ai, 1993) to understand the processes involved from the point of view 
of single systems. 

Apart from these more analytical applications of the quantum jump approach also its 
numerical usefulness has been demonstrated , for example, in numerical simulations of laser 
cooling experiments in two or three dimensions (see, for example, (Castin and M0lmer, 
1995; Marte et ai, 1993a; Marte et ai, 1993b)). It is the fact that the quantum jump 
approach allows the description of the system using a wave function instead of the density 
operator that made these investigations possible. Both the gain in computational speed 
and the saving in memory space is considerable as in a quantum jump simulation only A^ 
differential equations have to be propagated instead of A^^ in the density operator simulation. 

The quantum state diffusion model, apart from its importance in fundamental issues 
such as the measurement process or intrinsic decoherence, is now widely used to investigate 
the transition between classical and quantum behaviour and in the field of quantum chaos 
(see, for example, (Rigo and Gisin, 1996; Spiller and Ralph, 1994; Brun et al, 1996)). As 
we have already mentioned in Section IV the quantum state diffusion model also exhibits 
interesting localisation properties in both position and phase space (Gisin and Percival, 
1993a; Gisin and Percival, 1993b; Herkommer et al, 1996; Percival, 1994a) which can be 
used to implement very fast simulation algorithms (Schack et ai, 1995; Schack et ai, 1996a). 
Similar localisation properties exist also for variants of the quantum jump approach (Holland 
et al, 1996). 

A more exotic application of the quantum jump approach, or more precisely of the 
experiments in which quantum jumps were observed, is the fact that these experiments can 
provide a "perfect" random number generator. This is because these experiments allow the 
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observation of single quantum jumps which occur at absolutely random times due to the 
fundamental indeterminacy of quantum mechanics (Erber and Putterman, 1985; Erber et 
al, 1989). Whether this idea is useful is, however, doubtful although in principle there could 
be applications, e.g., in cryptography. 

F. The spectrum of resonance fluorescence and single system dynamics 

So far our examples of the quantum jump approach were limited to the investigation 
of the photon statistics of the radiation emitted by a system and of the internal dynamics 
of a system conditioned on the measurement record observing the radiation emitted by the 
system. As a further application of the quantum jump approach we would now like to 
consider the spectrum of resonance fluorescence of a three-level atom in the V-conflguration 
as shown in Fig. 1. This system, whose photon statistics we already discussed in the 
context of bright and dark periods exhibits interesting features in the spectrum of resonance 
fluorescence on the strong ^^ 1 transition. In the following we will discuss briefly the 
ensemble behaviour of the spectrum and then show how we can understand this behaviour 
from the point of view of the single system dynamics. We will closely follow the analysis 
in Hegerfeldt and Plenio (1995b) and Plenio (1994,1996). We consider the system shown in 
Fig. 1 and assume that the Rabi frequencies Qi of the lasers and the decay constants Ta of 
the two upper levels satisfy the following conditions. 

• The Rabi frequency of laser driving the ^^ 2 transition is weak, that is 

^^^i n, + {A,-A,y ^'''^ 

which simplifles for fii ^ Ai, A2 to 

^l « (0 . (223) 

• Spontaneous emission on the 2 ^-> transition should be negligible, that is 



These conditions have the following interpretation. If Eq. (222) is satisfied the V system 
exhibits long bright and dark periods as discussed in Eqs. (195) -(206). If Eq. (224) is satis- 
fied stimulated transitions from level 2 to level are much more frequent than spontaneous 
emissions on the same transition. Under the assumption of Eqs. (222)- (224) the spectrum 
of resonance fiuorescence takes on the following approximate analytical form 

1 /'''° 

S{A) = -^Re / dr{aio{r)a{0)),se''^^ 
npii Jo 

= Scoh{A) + SmoIIow{A) + Speak{A) (225) 

with 

B = 6Tl^-2Ql-2Al , (226) 

C = nt + 2nlTl^ + 9Tt^ + Al + 2Alnl-6Tl^Al , (227) 

D = Tl^{nl + 2Al + 2Tl^f , (228) 

^ _ , (A? + Tj,) [(Ai - A,)^ + r^,] mi - 4AI + 4Ai A,)^ + 16A^r^,] 
" f]2j^2r^^[^2^2A2 + 4A2-4AiA2 + 4rfi]' ' ^ ' 

2QlQlru{2Al + AAl-AAiA2 + Vtl + 4rl-^) 
" ~ mi - 4 Ai + 4Ai A2)2 + 16 Airf J {Qj + 2A2 + 2 T?,) ^^^°^ 



and 



''^''"^^"-^^>-2{A^ + BA^ + CA^ + D) ' ^^^^> 

nS,,ak{A)=^^^. (233) 

The contributions Eqs. (231) and (232) just represent the well known Mollow triplet and 
the Rayleigh peak (Mollow, 1975; Mollow, 1972a; Mollow, 1972b; Mollow, 1969; Mollow, 
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1981). These contributions are expected and are well-understood. The third contribution 
however, a narrow Lorentzian, is a new feature in the spectrum of resonance fluorescence. 
In Figs. 26 and 27 we see the spectra for the cases of strong driving of the ^^ 1 transition 
and for a medium strong laser on the same transition. One clearly observes the narrow 
peak , which should not be confused with the Rayleigh peak, in both spectra. In the 
following we will focus our attention to this new feature and use it as an example for the 
application of the quantum jump approach for spectral information and to illustrate how 
ensemble properties can be understood better from a single system point of view. One 
might, for example, be interested in the following question. What happens if we observe the 
spectrum of resonance fluorescence in a bright period exclusively? Experimentally this may 
be measured by triggering the spectrometer with the broadband counter. The spectrometer 
will be opened if we detect photons in the broadband counter at a sufficiently high rate and 
it will be closed when we fail to find photons for a certain time Tq. Using this time constant 
To we have a means either to observe the spectrum in a bright period (To ~ lO/Fn) or with 
no restriction on the emission times (Tq = oo), i.e., the ensemble spectrum. Changing Tq we 
can continuously switch between the two regimes. The spectrum of resonance fluorescence 
can be simulated using the quantum jump approach (Hegerfeldt and Plenio, 1996). In Fig. 
28 the effect of the change in To is clearly visible. For small Tq the amplitude of the narrow 
central peak becomes small. However only the central part of the spectrum is strongly 
affected while the wings of the spectrum are more or less independent of Tq. This is an 
example of conditional fluorescence spectra and it shows that spectra can be dependent on 
the conditions imposed on the photon statistics. Here we have shown the effect for the 
spectrum of resonance fluorescence but similar results may be obtained by investigating 
the absorption spectrum of a weak probe beam on the strongly driven <-^ 1 transition 
(Plenio, 1996). Again a narrow peak is observed which vanishes if we try to measure the 
same feature in a bright period of the system. After we have investigated the spectrum of 
resonance fluorescence of the V system and some of its properties, especially the occurrence 
of the narrow peak, it is now instructive to see how we can understand this feature from 
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the photon statistics of the single atomic system. The photon statistics are provided by 
the quantum jump approach. Exphcit results have been given in this section in Eqs. (195) 
-(206). 

Now that we know that the V configuration exhibit light and dark periods for those 
parameters for which the narrow peak in the spectrum appears we proceed with a somewhat 
simplified model of resonance fiuorescence of the V-system. We assume that the lengths of 
light and dark periods obey exactly Poissonian distributions, e.g., the probability density 
that a light period has the length t is 

h{t) = ^e-'/^^ (234) 

and for a dark period to have a length of t is 

lD{t) = ^e~'/^-. (235) 

These probability densities have been derived from a rate equation model of the time evolu- 
tion by Cook and Kimble (Cook and Kimble, 1985). Additionally we assume that in a bright 
period the system behaves exactly like a two-level system made up of the two levels and 
1. This assumption has to be checked and turns out to be good in the case of the spectrum 
of resonance fiuorescence treated here. However, it is not always correct as a careful anal- 
ysis of the absorption spectrum on the strong ^^ 1 transition proves and deviations from 
this idealised assumption may lead to significant contributions in the absorption spectrum 
(Plenio, 1996). 

It is well known (Loudon, 1983) that the stationary spectrum of resonance fiuorescence 
of such a two-level system is given by 

5(2)(A) = hm ^ f\ rdt2{E^-\h)E^+\t2))e~'^^'^~'''^ , (236) 

T^oo I Jo Jo 

where E'^~\t) and E^'^\t) denote the negative and positive frequency part of the electric 
field operator and C is chosen in such a way that with 

E^-^\t) r^ aoi{t) (237) 
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one obtains 

2r 



/ I /'OO 

5(2)(A) = ^iRe / dTe-^^^{a,,ir)aom)ss ■ (238) 

n Jo 

For the nonnormalized resonance fluorescence spectrum of tlie two-level system one then 

obtains 

c.(2). . ^ ^ ^n^l [1 Tnnl{nl + 2A' + 8Tl,) 2(Af + T^,) 1 

^^ fi2^2(A2 + r2,) \7r2(A6 + 5A4 + CA2 + D) ^fi2^2(A2 + r2,)'^^^/ ' 

(239) 

where the constants B, C and Z^ are given in Eqs. (226)- (228). As the light emitted by the 
atom switches on and off due to the light and dark periods we assume that the electric field 
radiated by the 3-level configurations is given by 

EW(t):=EW(t)/(t) , (240) 

where f{t) is a two state jump process with values and 1. The probability density for the 
length of a period where f{t) = is given by Eq. (234) and that for f{t) = 1 is given by 
Eq. (235). Therefore we have to substitute 

a,,{t):=a,j{t)f{t) (241) 

in Eq. (238) and expect the spectrum of the three-level configurations to be given by 

5(3)(A) = ^iRe / rfre-^-((aio(r)aoi(0))),, , (242) 

TT Jo 

where ((.)) denotes both the quantum mechanical average as well as the stochastic average 
over all realizations of the process f{t). This can be simplified to 

9 P POO 

^(3)(A) = ^Re / rfAe-^^^((7lo(r)(7ol(0)),,(/(r)/(0)),toc/^. 
vr JO 

/OO 
dA' S^^\A - A')k{A') (243) 

-OO 

with 



1 r'^ 
k{A) = - / dre-'^^{f{T)fm,toch. 
n JO 

Tl r Tr .,.. T, 



where 

Because both Tp and T/^ (Eqs. (204)-(206)) are assumed to be much longer than the mean 
emission time of a two-level system, which is of the order of Ff/, we can deduce from this 

._^{C(A, + ^.-,A,.^^^ii|±S^r,A)} ,,246, 

This expression has to be compared with the results Eqs. (225)-(233). In fact inserting Eqs. 
(204)- (206) into Eq. (246) yields an expression in very good agreement with the spectra of 
the V-system. The width Tp of the resulting narrow peak is 

r, = -^ + ;! . (247) 

The amplitude Ap of the narrow peak in the normalized spectrum is then given by 

Now the interpretation of the narrow peak is obvious. The stochastic modulation of the 
resonance fluorescence due to dark periods leads to a partial broadening of the Rayleigh 
peak. The small width of the additional peak in the resonance fluorescence spectrum is then 
understood from the fact that the correlation time Tc of the random telegraph process f(t) 
which simulates the light and dark periods is very large. In fact it is easy to show that 

1 1^-1 



T ' T / - (249) 

This results in an extremely narrow distribution in frequency space with a width Tp = t~^. 
It is this structure which is observable in the spectrum of resonance fluorescence. It is 
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interesting to note that the narrow peak is not easily interpreted in a dressed-states picture. 
Indeed, in secular approximation and for Ai = A2 = one obtains a zero weight for the 
narrow peak. Even if we tune the laser on the ^^ 2 transition to resonance with one of 
the dressed states, i.e.. A2 = ±f2i/2, the weight of the narrow peak comes out much too 
small. In fact it would then be predicted to be proportional to Q"^- The narrow peak found 
here has clearly a different origin than the line narrowing effects found by others (Manka et 
ai, 1993; Narducci et al, 1990b; Narducci et al, 1990a) where the systems do not exhibit 
bright and dark periods in their resonance fluorescence but simply a decrease in intensity. 
The existence of bright and dark periods in the resonance fluorescence leads to a narrow 
peak in the spectrum but it should be noted that the converse is not necessarily true. There 
are situations (eg a laser driven TLS in a squeezed vacuum (Swain, 1994)) in which the 
spectrum of resonance fluorescence exhibits narrow peaks but where the photon statistics 
does not show bright and dark periods. The reason for that can be found in the fact that the 
photon statistics is governed by the population decay rates while the spectrum of resonance 
fluorescence is derived from the g^^\t) correlation function which is strongly influenced by 
the decay rates of the coherences in the system. In the case of a TLS in squeezed light the 
coherences have a slowly decaying component while the population decay rate is still large. 

G. Spontaneous emission in quantum computing 

As a last application of the quantum jump approach we would like to investigate the 
influence that spontaneous emission has on the function of a quantum computer. Quantum 
computing is an idea that has attracted enormous interest in the last two years. It was 
elevated from the obscurity of theoretical idealisation to possible practical applications by 
the discovery of an algorithm by Shor (1994) (see also Ekert and Josza (1996) and references 
therein) that allows the factorization of large numbers in polynomial time on a quantum 
computer as compared to the exponential time required on a classical computer. However, 
this achievement in computing speed is only possible due to the massive use of the superposi- 
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tion principle in quantum mechanics. The basic idea is that a qubit (a two level system) can 
exist in a superposition of the two values and 1. N qubits can then exist in a superposition 
of 2^ values. These values can then be manipulated by a series of unitary transformations. 
A final readout can then provide us with information about global properties of the func- 
tion implemented by the unitary transformation. Such a global property of the function 
is, for example, its period which a quantum computer determines by performing a discrete 
Fourier transform, something which can also be implemented on a quantum computer (Cop- 
persmith, 1994). In the course of its time evolution (computation) the quantum computer 
evolves into a highly entangled state. However, it is known that any entangled state is 
very sensitive to dissipation. Therefore one expects that the quantum computer is highly 
sensitive to spontaneous emission and other sources of dissipation. This is in fact the case 
and currently research in quantum error correction methods has concentrated on attempts 
to find methods to correct for these errors (Shor, 1996; Steane, 1996; Calderbank and Shor, 
1996). The quantum jump approach is ideally suited For the investigation of this problem 
because it is able to describe single runs of a quantum computer rather than an ensemble of 
quantum computers as in the Bloch equation description. We will illustrate the problems 
caused by spontaneous emission in quantum computers by examining the example of the 
discrete Fourier transform mentioned above. There are two effects contributing to the deco- 
herence of the quantum computer. One is the obvious fact that a spontaneous emission will 
destroy at least part of the coherence in the quantum computer. The second decohering ef- 
fect, however, originates from the conditional time evolution between spontaneous emissions 
(Plenio et ai, 1996b). We have learnt above that this time evolution is actually different 
from the unit operation because even the non-detection of a photon represents a gain in our 
knowledge about the system. Therefore the wavefunction of the system, which represents 
our knowledge of the system, has to change. This leads to a distortion of the time evolution, 
which will then affect the result of our calculation. In Figs. 29 and 30 we simulate a quantum 
computer (we do not go into detail concerning its implementation here) that calculates the 
discrete Fourier transform of a function that is evaluated at 32 points (Plenio et ai, 1996b). 



The resulting square modulus of the wavefunction of the quantum computer is compared 
to the exact result obtained from a absolutely stable quantum computer. The function on 
which we perform the DFT is given for definiteness in this example by f{n) = Ss^(^nmodw) for 
n = 0, 1, . . . , 31. One can implement the Hamilton operators (in the Lamb-Dicke limit) for 
all the necessary quantum gates in a linear ion trap (Cirac and Zoller, 1995) to realize this 
DFT. In addition to the coherent time evolution, possible spontaneous emissions from the 
upper levels of the ions are taken into account but all other sources of loss are neglected. 

In Fig. 29 one emission has taken place during the calculation time of the quantum 
computer. If we compare the resulting wavefunction with the correct wavefunction, we 
observe a marked difference between the two. In Fig. 30 we show the wavefunction of 
an unstable quantum computer which has not suffered a spontaneous emission during the 
calculation of the DFT. We clearly see that even when no spontaneous emission has taken 
place, the wavefunction of the quantum computer differs substantially from the correct result. 
This difference becomes stronger and stronger the larger the ratio between the computation 
time T and the spontaneous lifetime Tgp of the quantum computer becomes. Therefore the 
wavefunction of the quantum computer will be sufficiently close to the correct result only if 
the whole computation is finished in a time T that is much shorter than the spontaneous 
lifetime Tgp of the quantum computer. 

The fact that even one spontaneous emission will usually make the result of the quantum 
computation completely incorrect can then be used to derive stringent upper limits to the 
numbers that can be factorized on a quantum computer (Plenio et al, 1996a; Plenio et ai, 
1996b). This again shows that knowledge of the single system behaviour gained from the 
quantum jump approach gives us useful new insights into important properties of quantum 
systems. 
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VI. CONCLUSIONS 

Recent work in quantum optics has forced us to re-examine the dynamics of individual 
quantum systems in which single realizations (single atoms or trapped ions, single cavity field 
modes and so on) are described in quantum mechanics. In these situations, the dynamics 
is always dissipative, and leaves a record of its history accessible in the wider world of the 
outside environment. If this record is read, so that we acquire specific information, then 
we can associate a specific quantum trajectory to that conditional record. In this way we 
"unravel" the dissipative master equation into a family of records. We have reviewed the 
new technique developed to describe this unravelling which go under the name of quantum 
jumps, Monte Carlo wavefunction simulations and so on. We have further demonstrated how 
they can be used to describe entirely non-classical behaviour in a wide range of situations 
in quantum optics. Future applications will surely emerge from these powerful approaches. 
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FIGURES 



1 




FIG. 1. The V system. Two upper levels 1 and 2 couple to a common ground state 0. The 
transition frequencies are assumed to be largely different so that each of the two lasers driving the 
system couples to only one of the transitions. The 1 ^^ transition is assumed to be strong while 
the 2 •i-^ transition is weak. 
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FIG. 2. Recorded resonance fluorescence signal exhibiting quantumm jumps from a laser excited 
'^'^Mg^ ion (Thompson, 1996). Periods of high photon count rate are interrupted by periods with 
negligible count rate ( except for a unavoidable dark count rate). 
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FIG. 3. Time-evolution of the population in the strongly-fluorescing level 1 of the three-level 
ion shown in Fig. (1). The lifetimes 7^^ and ^2 ^^^ marked on the time axis. What is crucial 
here is the "hump" Apu : this is a signature of the telegraphic nature of the fluorescence. 
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FIG. 4. A few periods of bright and dark sequences in the fluorescence intensity / from a 
three- level system. The bright periods last on average Tb, and the dark periods T^. 
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FIG. 5. Atom+field states used to describe the survival of a three-level system in an un-decayed 
state. The numbers Ug, riw specify photon numbers driving the strongly allowed transitions ^^ 1 
and the weak transition ^^ 2. Note that fluorescence takes the system out of the three atom+field 
states. 



109 



W(n,T) 



1/3 



Off 



An- 10 




On 



number of photocounts n 



FIG. 6. Macroscopic photocount fluctuations revealed in photo count distribution W{n, T) of 
counts of the allowed ^^ 1 transition fluorescence from a saturated V - configuration. One either 
registers a large number of counts in the time interval [0, T] (On) which is a sign of a bright period, 
or one registers no counts at all (Off) which shows that one is in a dark period. After Schenzle 
and Brewer (1986). 
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FIG. 7. Schematic picture of the heterodyne detection scheme. Cavity A emits a weak signal 
that we mix with the signal from the local oscillator cavity B. We measure the difference in the 
counts in the two counters that detect the photons that leave the beamsplitter. 



Ill 



/-A 



H. 



V-7 



H 



R 



£ 



Counter 



H 



B 



V-7 



FIG. 8. Schematic picture of the experimental situation envisaged in the theory of "cascaded" 
quantum systems. System B is driven by the quantum system A. The counter registers a super- 
position of both fields, one emitted from system A and the other emitted from system B. After 
Carmichael (1993b.) 
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FIG. 9. Ensemble averaged time evolution for the expectation value (0-3) = (|1)(1| — |0)(0|)/2 
(inversion in the two-level system). The dotted line shows a sample of 250000 trajectories obtained 
by the fourth-order Monte-Carlo method (17 = A, ^6t = 0.1 A^^ and zero detuning). It is hard 
to distinguish the dotted line from the solid line showing the analytical result. The dashed line 
shows a sample of 250000 trajectories obtained by the first-order Monte-Carlo method for the same 
parameters. From (Steinbach et al, 1995a). 
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FIG. 10. A schematic representation of a possible experimental setup for the measurement of 
conditional spectra. The spectrometer occupies a solid angle Qs while the broadband counter 
occupies Qb- The broadband counter performs frequent measurements while in the spectrometer 
only one measurement at a late time T is performed. 
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FIG. 11. The upper state population pn of an ensemble of two-level systems which is driven 
by a laser with Rabi frequency $7 = 5 F and vanishing detuning A = 0. The system exhibits some 
oscillations and then approaches a nonzero steady state value. 
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FIG. 12. The time evolution of the upper state population pn of a single driven two-level 
system. As in Fig. 11 the two-level system is driven by a laser with Rabi frequency fi = SF and 
vanishing detuning A = 0. The system starts a Rabi oscillation which is then interrupted by a 
quantum jump (detection of a photon). After the jump the system is reset in the ground state and 
a new Rabi oscillation starts. (After (Garraway et al., 1995)) 
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FIG. 13. The ensemble result of Fig. 11 is compared to the average over N = 100 and N = 1000 
realizations generated using the quantum jump approach. The fluctuations of the averages become 
smaller with increasing N and the ensemble average gets approximated more and more closely. 
The Rabifrequency is il = 5r and the detuning is A = 0. 
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FIG. 14. The next photon rate for different counter efficiencies /3 = 1,/? = 0.1 and /3 = .0049 
and the any photon probabihty which is obtained in the hmit (3 = 0. The Rabifrequency is fi = 5 T 
and the detuning is A = 0. We have normahzed the pfotted functions to have a maximum value 
of unity. 
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FIG. 15. The delay function Po{t) of the ^-system describing the probabihty that af- 
ter an emission at t = no other emission has taken place until t. Parameters are 
ill = 2rii,r22 = 0.35rii,Ai = A2 = and T22 = 0. One observes a slowly decaying tail of 
Po{t) indicating the possibility of dark periods. 
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FIG. 16. J, ,rp , representing the average rate at which quantum jumps will be observed in the 
three-level V-system, as a function of the detuning A2 on the weak transition. The parameters 
are Vti = 10rii,r22 = CSFujAi = and T22 = 0. One observes that the maximum quantum 
jump rate is achieved for detunings A2 = ib$7i/2 illustrating the Autler-Townes splitting of the 
ground state. 
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FIG. 17. The intensity correlation function ^'^'(r) of a y - system for the same parameters as 
in Fig. 15. One clearly observes that g^'^'{t) first falls off quickly to a value of around 1.2 for times 
around r ~ Srj"^ and then starts to fall off slowly towards the stationary value of 1. 
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FIG. 18. The time evolution of the population pn of the rapidly decaying level in the V system. 
Periods where the time evolution exhibits rapid Rabi oscillations interrupted by quantum jumps 
can suddenly stop and lead to periods of no Rabi oscillations and no jumps. The parameters are 
ni = 2Tu,n2 = 0.2rii,Ai = 0,A2 = andr22 = 0. 
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FIG. 19. The same parameters as in Fig. 18 but now for the population p22 of the metastable 
state 2. One clearly observes that during a dark period the population evolves smoothly into the 
metastable state 2. 
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FIG. 20. For the same parameters as in Fig. 18. The population of the unstable state 1 in 
a dark period. The population does not jump out of level 1 but evolves continuously. However, 
at the end of the dark period a jump occurs which is due to an emission for level 1 because we 
assumed the shelving level 2 to be stable. 
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FIG. 21. Single realizations of a driven two-level system whose resonance fluorescence is ob- 
served in homodyne detection. The Rabi frequency is fi = 4 Tn and the amplitude of the local 
oscillator is a) a = 0, b) a = 0.5, c) a = 1 and a = 10. With increasing a jumps become more 
frequent and smaller in amplitude. From (Granzow, 1996). 
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FIG. 22. A cavity with decay rate 2T is prepared in an initial state |^) = (|0) + |10))/2. We plot 
the expected photon number of the state of the cavity as a function of time. Before we observe the 
first photon outside the cavity, the expected photon number decreases. The first jump increases 
the expected photon number because we now know that the state has to be |9). Subsequently each 
photo detection decreases the expected photon number by one. 
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FIG. 23. Revivals in the inversion (a^) = (|1)(1| — |0)(0|)/2 of a two-level atom in a cavity with 
an initial field prepared in a coherent state \a) with a = 4. The parameters of the simulation were 
A = 0, g/j = 100, where g is the atom-field coupling constant. 7 is used to scale time while in 
the next figure it is the decay constant of the cavity. 
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FIG. 24. Revivals in the inversion (a^) = (|1)(1| — |0)(0|)/2 of a two-level atom in a cavity 
with an initial field prepared in a coherent state \a) with a = 4. Even a modest decay rate 
of the cavity leads to a rapid destruction of the revivals. The lower figure (b) shows the re- 
vival. It was obtained by a quantum jump simulation using 320000 runs and is indistinguishable 
from the numerical integration of the master equation. The parameters of the simulation were 
( A = 0, g/7 = 100, 7(5t = 5 X 10"^ .). 
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FIG. 25. For the same parameters as in Fig. 24 we plot a single realization of the time evolution. 
The diamonds mark the instants where a photon was detected outside the cavity. Note that the 
revivals in a single realization survive. The quantum jumps lead to a phase jump of the time 
evolution of the inversion of the atom which leads to a destruction of the revivals in the ensemble 
average. From (Garraway, Knight and Steinbach, 1995) 
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FIG. 26. The spectrum of resonance fluorescence on the ^^ 1 transition. The parameters 
are f^i = 6rii,r22 = 0.4 Fn and Ai = A2 = 0. One clearly observes the Mollow triplet and the 
additional narrow peak in the spectrum. From (Hegerfeldt and Plenio, 1995b). 
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FIG. 27. The spectrum of resonance fluorescence on the ^^ 1 transition. The parameters 
are ili = 2rii,il2 = 0.2 Fn and Ai = A2 = 0. Due to the weaker driving as compared to 26 
the sidebands in the Mohow triplet are no longer resolved. The additional narrow peak is again 
observable and has a higher relative weight. From (Hegerfeldt and Plenio, 1995b). 
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FIG. 28. Simulation of the conditional spectrum of resonance fluorescence for Tq = oo and 
To = 50 rj"^ . One observes that only the narrow central peak is affected by the conditioning of the 
photon statistics while the wings of the spectrum are essentially independent of the choice of Tq. 
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FIG. 29. Results of a discrete Fourier transform (DFT) of a function f{n) = Sgtnmodw) with 
n = 0,1,... ,31. The solid line is the result for a quantum computer with stable qubits and 
represents the correct result. The dashed line shows the result of the same computation using a 
quantum computer with unstable qubits, one of which has suffered a spontaneous emission during 
the calculation. The results clearly differ and show the impact of a single spontaneous emission on 
a quantum computation. For the parameters chosen on average the quantum computer will suffer 
one emission per DFT, ie Tgp = T in this case. From (Plenio et al, 1996b). 
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FIG. 30. The same quantum computation as in Fig. 5. The sohd hne again represents the 
result using a quantum computer with stable qubits, while the dashed line shows the result using a 
quantum computer with unstable qubits. This time, however, the unstable quantum computer does 
not suffer an emission during the whole calculation. Again the results differ illustrating the impact 
of the conditional time evolution between spontaneous emissions. From (Plenio et al, 1996b). 
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